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ABSTRACT 


This thesis describes a study of the rotationally invariant 
part of the electron-coupled nuclear spin-spin interaction, other- 
wise known as the Ruderman-Kittel interaction, in 8 (or white) tin. 

It is a continuation of the NMR study initiated in our laboratory 
by Smith (1972), in which the anisotropic nuclear spin interactions 
in 8-tin are averaged by rotating a sample of tin powder at high 
speed using an air turbine. Provided the axis of rotation makes the 
so-called magic angle (eosmeiny3) with the external magnetic field, 
the broadening associated with the anisotropic interactions is com- 
pletely removed and the NMR lineshape is governed by the Ruderman- 
Kittel interaction together with spin-lattice relaxation. The 
theoretical lineshape is synthesized and the strength of the Ruder- 
man-Kittel interaction is obtained by fitting the theoretical line- 
shape to the experimental resonance. 

Improvements made to the air turbine designed by Smith enabled 
us to record the Sn!!7 NMR signal with a greatly improved signal-to- 
noise ratio. The theoretical lineshape calculations made by Smith 
have been extensively revised and fewer approximations have been made. 
Also, a new fitting procedure has been developed which enables a 
theoretical lineshape which is computed numerically to be fitted to 
an experimental rneetene: 

Agreement between experiment and a theoretical model, which as- 


assumes that only nearest neighbour spin-spin interactions are signi- 
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ficant, is very poor. Similarly, a model in which only second near- 
est neighbour interactions are non-zero, has been shown to be unten- 
able. However, very good agreement has been obtained with the so- 
called Ruderman-Kittel model and we find that the strength of the in- 
teraction between nearest neighbour Sn!!7 and Sn!!9 nuclei is 3.04 + 
0.3 kHz. A comparison is made with other published values of this 


quantity. 
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I INTRODUCTION 


I(i) Generai 


It has been known since the early work of Bloembergen 
and Rowland (1953) that the NMR linewidth of nuclei with I=1/2 
cannot be wholly accounted for by direct nuclear magnetic dipole- 
dipole interactions and by a lifetime-limiting effect associat- 
ed with spin-lattice relaxation. By measuring the T1293 and 
T1295 linewidths in various samples of thallium metal with dif- 
ferent isotopic composition, Bloembergen and Rowland (1955) showed 
that an indirect interaction occurs between nuclear spins via 
the intermediary of the. conduction electrons. In particular, 
they showed that the indirect interaction between two spins a 
and a consists of a scalar or rotationally invariant part of the 
—— fone = given the name of pseudo-exchange inter- 
action, together with a smaller anisotropic part, sometimes cal- 
led the pseudo-dipolar interaction. A theory of the rotation- 
ally invariant part has been given by Ruderman and Kittel (T954)% 
and indeed this interaction is also often called the Ruderman-Kit- 
tel interaction, a name which we shall adopt in this thesis. 

The Ruderman-Kittel interaction between unlike neigh- 
bours contributes to the second moment of the NMR line (Ruderman 
and Kittel, 1954) and hence broadens it. On the other hand, the 
Ruderman-Kittel interaction between like neighbours leaves the 
second moment unaffected but increases the fourth and higher (even) 


moments. This implies a sharpening of the central part of the 
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resonance, a phenomenon called exchange narrowing. 

The study of the indirect nuclear spin-spin interaction 

is important on two counts: 

1) It provides an understanding of the factors which contribute 

to the width and shape of NMR lines in general. 

2) It provides information regarding the hyperfine interaction be- 
tween the nuclear spins and conduction electrons in metals. 

There is some uncertainty in the magnitude of the Ruder- 
man-Kittel interaction in the case of white or 8-tin, since all 
the methods used so far have required a knowledge of the second moment 
of the tin resonance. Values for the magnitude of the Ruderman- 
Kittel interaction between Sn!!7 and Sn!!9 nuclei which occupy 
nearest neighbour sites in 8-tin range from 2.0 + 0.5 kHz obtained 
by McLachlan (1968) to 4.1 + .3 kHz obtained by Alloul and Deltour 
(1969). A method has recently been described by Smith (1972) and © 
uses the so-called magic angle, high-speed rotation technique. 

This involves rotating the sample at high speed Gerace at 
several kHz) about an axis making an angle cosa (473) with the 
external magnetic field. 

In this thesis we describe some additional work done on 
this problem. In particular, substantial refinements to the theory 
given by Smith (1972) are presented together with a description 
of improved rotors which enable a much better signal-to-noise ratio 
to be achieved. Also, a better procedure for fitting theoretical 
lineshapes to experimental resonances is described- 


In sections I{ii) and I(iii) of this thesis, a simple pic- 
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ture of NMR and a discussion of the magic angle, high-speed rota- 
tion technique are given. Section II contains the theoretical 
background for this thesis. In section III a discussion of the e- 
quipment and experimental method are given. In section IV the fit- 
ting procedure is described and the results are presented and dis- 


cussed. 
I(ii) A Simple Picture of NMR 


Purcell, Torrey and Pound, (1946) and Bloch, Hansen and Pack- 
ard (1946) carried out the first spin resonance experiments on nuclei 
in liquids and solids. When a nucleus with spin quantum number I 
is subjected to a steady magnetic field Ho? there are, 21+) equal-— 
ly spaced magnetic energy levels, with a separation 

AE = WH/I ~ yhH, (1) 
between adjacent energy levels. Here u = yhI is the maximum z 
component of the nuclear magnetic moment (2 is taken to be in the 
same direction as H) and y is the gyromagnetic ratio of the nu- 
cleus. If a group of nuclei which interact negligibly with each 
other are allowed to come to equilibrium with a heat reservoir at 
temperature T, a net magnetization, M in the z direction is es- 
tablished which obeys the well known Curie iaw M « 1/T. When 
rf-wmagnetic field of frequency 

w = yH (2) 

is applied perpendicular to Hs there is a net absorption of ener- 


gy from the rf field with a new equilibrium value of Me Because 
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of the finite lifetime of the magnetic energy levels due to the in- 
teraction with the heat reservoir, there is also absorption of en- 
ergy from the rf field for frequencies close to Wo For a group of 
spin 1/2 nuclei which interact only with a heat reservoir, the shape 
of the absorption spectrum is Lorentzian, the half-width of the spec- 
trum being equal to 2W/1, where W = eae is the transition rate 

and - is the spin-lattice relaxation. time. 

For solids, the direct magnetic dipole-dipole interaction 
between nuclei is important (Van Vleck, 1948). Classically, the 
dipole-dipole interaction can be considered in the ee way. 
Each nuclear magnet finds itself not only in the applied steady 
magnetic field Hee but alison assmal? local field Hoe produced 
by neighbouring nuclear moments. The magnitude and direction of 
Ayo differs from nuclear site to nuclear site, depending on the 
relative disposition of their magnetic quantum number m (where m can 
takeson Walues.-l, -l+b,.. 2. 1)2. The: distribution (0b Loe has a 
mean value of zero and a width of the order fone where r is the 
nearest neighbour distance. As the resonance frequency of each 
nucleus takes on different values due to the variation in aoe? the 
resonance will be further broadened (in addition to lifetime broad- 
ening). 

In a metal, another important interaction is that between 
nuclear spins and the electronic spins. A simple way to examine this 
is to consider the magnetic field produced by electronic spins. 


This rapidly fluctuating field can be divided into its time-average 


and the difference between the time-average and the value of the 
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field at a given time. The time-averaged field, which to a good ap- 
proximation is parallel to Hoe causes a shift AH of the resonance 
which is proportional to H, (and usually down field). This shift was 
first observed by Knight (1949) and the so-called Knight Shift K is 
defined as AH/H, - When the time-averaged distribution of the elec- 
trons about a nuclear site has lower than cubic (or tetrahedral) sym- 
metry, the Knight Shift is anisotropic, that is it depends on the 
orientation of the metal with respect to H,- It can be shown that 
the Knight Shift is a second rank tensor. This is a further source 
of broadening when the sample is in the form of a powder since the 
various grains of the powder will have different orientations with 
Tespectto Ho: The time-dependent magnetic field caused by the 
electronic spins gives rise to a spin-lattice relaxation. 

A second effect of the time-averaged magnetic field is the 
indirect interaction which was previously mentioned in section Wet). 
The existence of such an interaction can be visualized by the fol- 
lowing picture (Slichter, 1963). Consider two magnetic nuclei a- 
mongst non-magnetic nuclei. The effect of the spin of one of the 
nuclei is, in general, to make the nuclear site more favourable for 
an electron with a spin parallel to the magnetic moment of the nu- 
cleus. The zero states of the electrons are Bloch states and in 
order to increase the probability of an electron with a negative 
spin orientation being in the vicinity of a nuclear site of negative 
magnetic moment, it is necessary to mix Bloch states of negative spin 
orientation and with wavevector k greater than the Fermi wavevector 


k (because of the exclusion principle). These Bloch states are in 
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phase at the nuclear site but become progressively more out of phase 
with increasing distance from the magnetic nucleus. This implies 
that the conduction electron density oscillates about the average 
value. The magnitude of the oscillations decreases with increasing 
distance. Consequently, the second magnetic nucleus experiences 
a different time-averaged magnetic field, thereby giving rise to an 
interaction. 

There is also a nuclear quadrupole interaction but this in- 
teraction does not occur in tin as the magnetic isotopes of tin have 


a spin of 1/2. As a result we will not describe this interaction. 


I(iii) Macroscopic Rotation 


The first studies using the sample rotation technique were 
made by Andrew et al (1958) to confirm the theoretical predictions 
of Anderson (1954) and Pake (1956) that hindered rotation in a solid 
should not cause a reduction in the second moment of the NMR spec- 
trum. When they rotated the sample as a whole, Andrew et al found 
that the central part of the spectrum was narrowed but sidebands 
appeared in such a way that the second moment remained invariant. 
Narrowing by macroscopic rotation at the magic angle has been studied 
in more detail by Andrew and Newing (1958), Andrew and Jenks (1962) 
and Schwind (1967). It has been shown by Kessemeier (1967) that if 
the rate of rotation is much greater than the NMR linewidth of the 
stationary sample, the central part of the resonance is resolved 


from the sidebands, and only the time-averaged interactions need be 
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included in the Hamiltonian to study the central part. A comp- 
rehensive study has been undertaken by Andrew and Farnell (1968) in 
connection with the averaging of anisotropic interactions by rota- 
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LI? “THEORY 
II(i) General 


The Hamiltonian of a group of nuclear spins interacting with a 
"thermal bath' of conduction electrons in a eed can be written as 
/E = hE + HF + iQ (1) 
where fi—E involves only the nuclei, AF involves only the thermal 
bath and fiQ is the interaction between the nuclei and thermal bath. 


Q can be written as 


Kerr ert (2) 


=f 
Q ful hd q aa’ fufu 


/ 
where a,a@ are eigenvalues of E, and f,f are eigenvalues of F of 
degeneracy u and u respectively; Ka is an observable of the nuclei 


and Ht is an observable of the thermal bath. For a metal we have 


AQ = al.-2S,6(r,,-R.) + tb.1.- 2D, Sy + Ze, UsEL(r, ,)/r, 55) 
Bical) a als i imi Celie ala d iL 
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(re eRe) Cr Seo es 44 tand y. and y, are the gyromagnetic ra- 
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ttoseor thes ith nucleus. and ‘an; electron respectively. L(@®; 9) is 
the angular momentum operator for the &£th electron with respect to 
tnesith nucleus. as center tand Lig is the position vectonm for the 2th 
electron with respect to the ith nucleus (Das and Mahanti, 1968). 
The first term in (3) is the so-called Fermi contact interaction and 
the second term describes the dipolar interaction between the nuc- 
lear spins and the electron spins. 


The density matrix operator (Davydov, 1966 p. 42) for the 
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nuclei and the bath is p and is assumed to be of the form (Redfield, 


1957) 
Pea i Bee een (4) 
with 
Pel) a= exp (-RE/KT) /2, pxp (E/E), (5) 
u 
Redfield found that 
do/dt = (o(t+At)-o(t))/At = ifo,E+M+N] +R (6) 


Te At >> 7 SY atid M,NR << 1/1 < ki/h, where t can be thought 

& e S 
of as the classical correlation time for the random magnetic fields 
produced at the nuclear sites by the electrons. Redfield showed that 


the matrix elements of the operators M and N are 


Maat Ka fag” fufu 3 PCE) (7) 
Pee eee 
phe ay veccy™ 2 MN Gee ee i) ; (8). 
Also he found that 
MHD eae BR (9) 
and 
Ra = DK KG og (o-B) +5, 4-8) - (10) 
5g SK i aa = 5gik ok 5 glene epee 
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tage x sas (2 PCV e ued Hae I Gan rad 


and n,,(f) is the density of states for eigenvalue f and degeneracy 


u. In writing (10), Redfield assumed that only secular terms con- 


ie 7 


ne 


a. o 


% 7 D 


at: qq, v 


nae 1) 


oa 


yak 
7) 


y oad mos 


s sahE(t)e6€ whey +} 


Fou db 
el T\ IH 
i 
t a 
tS wi .® 
r Ta} 
. £5 
! 


, bp ' os te A 
a | 


r - 7 Ny 
‘ ted it oad 


s+ 


a Pee 


> 
" a _, 


Pot th 


a 4 
tributed to Maye that is only terms where a-a-B+8 = 0. However, 
for non-secular terms, Redfield showed that one must multiply the 
right hand side of (10) by a term of the form 

: , f. 
et A / 4 
Anofgg = ter A PFFIAE 137 (aa pe sat, (11) 
M can be identified with the Knight Shift operator (Kt 

in (2) is taken to be one of the components of the spin operator 
for one of the nuclei because of (3)) and N is the indirect spin- 

Cee ; ; : : / 
spin interaction operator (ignoring terms in (8) where q and q re- 
fer to the same nuclear spin). The R operator is the relaxation 
operator. 


Phoesxoperator: N(E+MtN) ‘cans therefore be expressed, as 


an yed: pts Oneal Zeal deat De gy) ole (12) 
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WieresweS the Knight shitt tensor, and Ji; and De areurespectively 


the indirect tensor and the dipolar interaction tensor between nu- 
cleist and. jo) Itotollows that (6);-can be writtem an ithe tor 


do/dt = i {o, (-2y,51;@Q*W-H, = Ti (L,5+2i,) 1) +Ro, (13) 
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11(@@i) “Effect of Macroscopic Rotation on K, J, and D 
Because of the short correlation time of the conduction elec- 
-11 
trons (t,. = 10 seconds) and the fact that the period of the macro- 


scopic rotation is typically me seconds, we can allow K, J and D in 


(13) to be time dependent. However, as pointed Oug “ing £0147) tire 
shape of the central part of the resonance is governed by the time- 
average of K, D and J. The R operator in (13) is-assumed to be un- 
affected by macroscopic rotation, since McLachlan (1968) found that 
the Ty of Sn/!9 in g-tin was essentially naomeriene of crystat 
orientation (the connection between Ty and R is established in sec- 
tion IT(iv))~ 

To find the time-averaged operators, it is noted that the in- 
teractions involving K, J and D are of the form A-I'B where A and B 
are time independent vectors. Because of the symmetry of the crystal 
Structure or B-tin, | must be symmetric matrix (Smith, 1972, p. 94). 
Following Andrew and Farnell (1968), we can decompose the tensor into 
two parts, T and T, where T is equal "to. (17.3) Tr(L) 1s and T is 2 
Peacevess Censor equal to 1)"- Zs Because i is symmetric 


*& * * 


where Cand Con are the direction cosines between the axes x, 
an 
y and z fixed with respect to the laboratory (represented in (14) 


by the indices a and 8) and the principal axes of the tensor. The 


* 
T,,'S are the principal values of the tensor T . Using the relations 


K, = AS + iAy and A_ = Ay - iAy plus similar relations for B, and B_, 


we obtain 
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We now consider. the. truncated part) .of (15), that is, we consider 
the terms Dp and E,. We can ignore the remaining terms when A and 
B represent: spin vectors since these terms give rise to weak reso- 


nances at frequencies 0, 2w,, and 3w,, resonances which we are not 


oO? 
interested in. We have 


trunc 
( 7g n°on’BnAaBg) 


2 Dene? 
BT, {CznAzB,+ (1/4) (A Be ty A0B,) (Cato 


2 | | 2 
BT (Cp, B+ 1/2) (A,B, TANB)) CTC ea 
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ET, (-1/2)C 57 (A,B, AUB 2A,B.), (16) 


z pea WAS awe 


In appendix I we show that the time-averaged value of Co 


oe 


* 


sesh ooeh : 2 
the. M2et Ce angle condition 1s\satisficd.) Since 7 wis traceless and Ca 


is 1/3, the right hand side of (16) reduces to zero. Thus we have 
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HJ 


i . (17) 


If truncation is justified, then the time-averaged interaction (when 


the magic angle condition is satisfied) reduces to 


AZ 


roa 
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be) 4 2 OD (18) 

For the dipolar interaction, A and B are spin vectors, and 
truncation is justified. Moreover, the dipolar interaction tensor 
is traceless, As a result, the time-averaged interaction ihe equal 
tol zero. 

The indirect interaction between nuclei can also be written 
in@the form AL By wheres Avand™B aremspinfoperators. OAs ay result 
the time-averaged indirect interaction becomes 


iG Co SOS Ok Alia UNG (19) 

isja Gla i<jecvion 

where ” = Ce eda elb 

Finally, the Zeeman interaction can be written in the form 

ACER However! AY is#arnuclear spins vector and) Bwis! the steadys mag- 

netic field Hy inwthisMeasere Nher stn of Hy and the magnetic: field 

caused by theselectronic™ spins is, to a‘eood approximation, -in the 

z direction (Abragam, 1961) and the interaction can be written as 
-rh kai i 20 
i Oe a) Ze Ogee ea) 

It follows that terms Ee through Le Ie CloanewzeLorile Onis case, 

and equation (17) obtained for the truncated interaction can there- 


fore be applied. Thus the time-averaged Zeeman interaction is 


~Dhiy; (1#K)T; 7H, (21) 


where Ko =.:(1/3)Tr(k) = (1/3) (K,+K,+ Kz). For $-tin K, = Ky and 


5 
K, = Ky = Ky because of the 4 axis of symmetry, 


In summary, the time-averaged value of A(E+M+N) is 


fi(E+M+N) = : hy; (A+K)1, Ho + J ear es (22) 


when the magic angle condition is satisfied. 
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II(iii) The Ruderman-Kittel Interaction 


The usual method of finding a theoretical expression for the 
Ruderman-Kittel interaction is to consider fiQ as a perturbation of 
A(E+F) in (1) and to find the second-order correction to the energy 
levels. The first order correction is responsible for the Knight 
Shift and will not be discussed here. The second order correction is 


id Ge 
AE = I BE /(BtE,-E Ey) (23) 


where | no’) is a product of a many-electron state In) Olmenern2 vi Ey? 
and a many-nucleus state |) of energy Ey? lie unperturbed state of 
the electrons and nuclei is | Oa). Examination of equation (3) shows 


‘> that is,it can be written as a sum 


that nO Ws ofethe form s(H ”) 
i en 
of terms, each one involving only one nucleus. By expanding fQ in 


terms of (H one obtains terms. of the form 


en) i 


| (na| (H | 0) |7/ (E+E, -E,-E,) (24) 


aes 


which involve only one nucleus. Since we are only interested in in- 
teractions involving pairs of nuclei such terms as (24) are omitted 
and we retain only the cross terms of (23), and thus 


AE =, FE (0a[ (Hyy)q [ne) (na | (Hoy); 100) / EgtEy-En-Ey). (25) 
JS cate) 


Because the electronic energy states are essentially continu- 
ous and the energy difference between the unperturbed and perturbed 
electronic states is usually greater than IEE yl» the denominator of 
(25) is approximated as E,-E,. Also, the interaction (How takes the 


form I*G., where G, does not involve the nuclear spin operator. Thus 
ed ecec sa 
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We could have obtained AE as a first order perturbation 
0a 


contribution if the term & had been added to (1) where 
cee) ee lee (eee) ork ts Zi 
H eff i<j ese ig) BIS tas oe 
and 


el esti 2{(0]G, ,[n) 16, ¢|0) + c.c.}/(E-E ). (28) 


H eee is equivalent to operator N of (6) if the terms of N which in- 
volve only one nucleus are excluded. We shall henceforth assume that 
newindirect. interaction as given bys (27). 


The right hand side of (28) can be written as (Mahanti and 


Das, 1968) 
(55) a6 jee + pe * Oey = a5) pe" oe 
where 
eee, = rB,((0|65—" In) @lG5qr 10) te.c.} 
eae = zB nf (O|6se"* [n) (nlo5a? 0) + Colcy,Pin) (n|csa” [0)} 
(J, a = IB fcolege’ In ) (nes? [0) siete} 
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Some of the cross terms involving the contact, dipolar and orbital 
contributions have been excluded from the right hand side of (29) 
since they are equal to zero when the orbital momentum is quenched 
(Mahanti and Das, 1968). 


If the electron energy surfaces are spherical, the term 


(2) 


(S55) 88 is dipolar in form and contains no diagonal components (Ma- 


(3 ? 


hanti and Das, 1968). The term (J, 58 » again assuming spherical 


(4) 


energy surfaces, has only diagonal components, while the term (J. 5 8p 


is composed of two parts, a dipolar part and a diagonal part (Ma- 


hanti and Das, 1968). These authors showed that for Rb and Cs, the 
&) and (J y 4) 


Concriputions of (J...) 7 tej}. LOethe diagonal terms or: Jie. 
Lee 1388 ; els) 
are small compared to the contributions of Gaae 
For Os ae we have 
area = aja; zt 2 Bt (nks|S, 6(x-R. 5) Inks) Gk 4 S(x-R. Bic) + 
BB nks nk 


een (30) 
where the sum IM SU) aS Over: OCCUpLed sbLOCNmeLeCtronestaces <1 Al 
nks 
SO, ni represents the band of the Bloch state in the first Brillouin 


zone, k is the wavevector and s is the electron spin of the electron 


State. The sum £,is over the unoccupied electron states and B is 


ak's’ 
equal to (E_ -E..-,). This expression has been rewritten as 
ks nks 
k-k 
GV, ee” oft gmk |§ (x) (nik) Gtk] 6 (x) [nk)e- iUCK)«R AH(SiS 4 \s)B(s|S, ‘|s) 


by various authors (see for example Slichter, 1963). Here Bae = R, - 
R_ and Inks) = |s) |nk). The probability of a state Inks) being oc- 
ae ia i ae 

cupied is f(nks) and of a state [nks) being unoccupied is 1-f(nks), 


where f is the Fermi-Dirac distribution function. It follows that 
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where the sums 2 and ~ are now over all states. Because f(nks) 

nks nks 
is approximately the same for spin up or spin down states, and matrix 
elements vary slowly with energy, we can replace f(nks) by f(mk) and B 
by (Eqym En! y+ (Slichter, 1963). As a result of this approximation, 


and since 


Bee ls a l=) e182) = TENS Ee = Tr(S,S4) = 64/2, (32) 


we have 


(1) : (1) 
COR up yma) rl aes 


= -(1/v)2 2 ET (kK E(ak) (1-£(nK) te" SO 
nk nk’ nn ee 7 pie 


J/(E-E ws, 33 
on fe (33) 
where I_,(k k) = -a.a ap (0) |*|u-¥(0) | and ~,(r) is the wavefunc- 

nn —’— i jz nk nk nk = 

Laon Of |nk) imvthescoordimate representarion. = Difrerent, authors have 
made various approximations in applying (33). In particular Ruderman 
and Kittel (1954) assumed that the electron energy surfaces are spheri- 
cal, and that 1 (kk) can be approximated by its value at the Fermi 
surface. Ruderman and Kittel found that 

J.. = (2/9n)y2y.y:h2m |p (0) | F (2k fre) (34) 

iy YeYiYj Ke sees 


4 * 
wmnere F(x) = ixcos(x) - Sin(x)!/x andm is the effective electron 


mass, The form of F(x) is shown in figure 1. A more general ex- 
pression was obtained by Mahanti and Das (1968), who lifted the re- 


striction that J i(k, K) be evaluated at the Fermi surface and that 
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Figure 1 Graphical Representation of F(x) 


The number n shows the position of x = 2k Riis where Ri; is the 


distance between nth nearest neighbours. 
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the bands be parabolic. They obtained 
UGE 1s 
gh) s fe/2(2n)°R, 7} Sdk m(k)k¢(nk,nk) sin(2kR:; .) (35) 
ij J 0 a ay | 


, 2 2 
where c = -(16/3)*17y y_y 2n4 and (nk nk) = lw (0) | |w+-(0) | 
5 i les i i nk nk 
Finally Roth, Zeiger and Kaplan (1966) lifted the restric- 
tion of spherical energy surfaces and derived more general expres- 


(1) (1) 
J 


sions for Je. 4 in particular theyetound that J; was proportional 


J 
-1 
to Se for directions. perpendiacular/to two) flat repions; of the Fer= 


i (1) =2 
mi surface and Ee was proportional to i for directions perpen- 


dicular to the axis of a ‘cylindrical region of the Fermi surface. 
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II(iv) Theoretical Lineshapes of Configurations of One, Two 


and Three Magnetic Nuclei 


In this thesis we shall be comparing the experimental Sn!!7 
resonance in a sample of B8-tin which is rapidly rotated at the ma- 
gic angle, with a theoretical lineshape in which only spin-lattice 
and Ruderman-Kittel interactions affect the lineshape. Because of 
the small abundance of magnetic isotopes in natural tin (0.35%, 
7.67% and 8.68% for Sn!15, Sn!!7 and sn!19 respectively (NMR Tables, 
Varian Associates, 5th edition, 1965)), we synthesize the theore- 
tical lineshape by combining, in ae proportions, the line- 
shapes, of nuclei which are interacting with zero, one or two other 
magnetic nuclei, ine probability of finding three or more magnetic 
RucCreI wim thegvicinity Olea mparticudar nucleus, asusmad |. 

Before proceeding, we shall discuss the application of equa- 
tion (6) to the particular problem of calculating lineshapes in §-tin. 
This equation was obtained by assuming that the terms M, N and R are 


much less than 1/t.. For §-tin the Knight Shift is approximately 


so 


s156 Womith,, 1972)) wath avresule that, for .acfield of approximately 

5 kG, the order of magnitude of M for 8-tin is 50 kHz. The magni- 
tude of the Ruderman-Kittel interaction for nearest neighbours is 

a few kHz. Finally, the magnitude of ae of the order of 6 kHz since 
T, is approximately 170 microseconds for Sn!!7(it will be shown in 
this section that R is Of the order tof 1/T). Since 1/t, is approx- 


a. vik y 
imately 10 seconds , M, N and R are much less than Ice 
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When non-secular terms are included in equation (6), the ef- 


fect of the terms Rs y for which (a-a+B-8) At > 1 can be ignored 


aBB i 
compared to the effect of terms R fad for which 
aa 
(a-o+8-B)At << 1 (36) 


because of the terms A ee {see equation (11)). As we neglect terms 
ote 

R dad Which do no& satisfy (36) when calculating the lineshapes, we 

Te) 
can replace A | u7-by. the; value 1 (the magnitude of A ,. ¢ approaches 

aaBB aaBB 
1 as (a-a+8-8)At approaches zero) and regain (6). 
The equation can be simplified by considering the high temper- 

ature limit and short correlation time approximation, that is when 

|E -E,|/h << wae kT/h, where EY and EY represent “eigenvalues of 

T Ai 

fi(E+M+N). We can then replace Ro = R(o-6 ) by R(o-o ) where R is de- 

: : : hw/2kT | : 

Tivedstrom Rk aby, using) ko (wissse ji, (Oe inv place-Otwy 17.0), ehe 
qq qq qq 
result being that Ridee = Re doc’ Also we can set ag = K 460) 
(Redfield, 1957). These approximations can be made since at room 

13 
temperature (the temperature of the 8-tin sample), kT/h = 4x10 ©sec- 
-1 -11 8 -1 
pnds: ©, ie = 10 seconds and JE Ey | /h = 10 seconds . We make two 
assumptions about Sess simplify thescalculations. The first as 
that k (0) = 0 unless q = aq’. Classically this means that. there is no 
qq 

correlation between the random magnetic fields at different nuclear 
sites and between the various components of the random field at a given 
nuclear site. This turns out to be a good approximation for g-tin (Win- 
ter, 1972). The second assumption is that Sa = k’ where k is inde- 


pendent of q. This is a good approximation since McLachlan (1968) 


found that qT) for Sn!!9 was essentially orientation independent. 
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Next we show that the resonance lineshape is unaffected 
by the operator N which appears in equation (6) provided N com- 
mutes with M. » the x-component of the magnetization operator. The 


expectation value of M,, is given by 


<M,> = Z (a [OM |) (37) 
and thus <M,> = d<M,> = E(aloM,Ja) (38) 
OQ 
dt | 


Substituting (6) into (38) we find that 


<M,> = (a [o ,E+M+N}M, |o) (39) 
As a result, the value of <M,> is not affected by N if Z(a|{o,NJM, |a) 
a. 
is zero. This is easily shown to be the case if (Msn = 0. We 
have 
Z (| oNM, -NoM, |) = Z (a |oM,N-NoM, |e) = Tr(oM,N)-Tr(NoM,) = 0 
Since the trace is invariant under cyclic permutation of o, M, and 
N. 
Finally we include the rf field (linearly polarized in the x 
direction). gif H : rsethe at fieldwmeteeract r0njuweshavemioreequation. (6) 
r 
5 ee i(E/-E jo Jet Rue 7 ‘(0 /-0 1 a 
ao a ar" ad "ga Cfo oy Sie aisis sys: 
(40) 
: “ 4 4, 
a Balt (OH, gle) Colt elo oy ; 
Here E, is an eigenvalue of the eigenstate |) for. the “operator 
A (E+M+N), and the subscripts refer to the different eigenstates of 
this operator. From (10) we find 
Z 
Rees (1/2h ) (25 pa fg Sy Qed 42-6 440 ) (41) 
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2 Dad, , / 

pk Col 15/8) (BIT; q 10) (42) 
The index j is summed over the nuclei and q takes the values x, y 
and z. In obtaining (42), we have made use of the approximations 
K) 
oie aes 


In the following subsections A), B), C) and D) we consider the 


involving Kad aq’ = 6 


spectrum of nuclei which interact with no unlike nuclei, one unlike 
and no like nuclei, two unlike and no like‘nuclei, and one like and 
one unlike nuclei respectively. For these cases we have Oey CtKH 


and Wa is the relaxation: rate of the Nth nucleus. 


A) No Unlike Nuclei Affecting the Resonance 


The lineshape function in this case has been shown previous- 
ly (Smith, 1972) to be given.by 
x(w) = Ariat) Wea (43) 
where Zy = ease a. We note that eauation (43) applies whatever 
the number of like nuciei which are interacting since eee where 


i and j are any two nuclei of a group of like nuclei, commutes with 


My = -Zy pfily,. the x-component of the magnetization operator for the 
Q 
B) One Unlike Nucleus (and No Like Nuclei) Affecting the Resonance 
For this so-called XA case we generalize the calculation pre- 


viously made by Smith (1972); we allow for different spin-lattice 


relaxation rates ah and ths for the X and A nuclei. 


AU 


| ak 4 aa | 
a i — hed 


Bis? * Sea 


ae ee ad Barres nen Wi 


ae 


‘tot taatxongea ol : 
~ 


Ge 
w , Pet) re 
ea¢ tabienco sw Or tone i a rh ebm Tai tire 


,2 { 
itiny ono , deboun sais oh at tw 


y 
Pape) 


The four eigenstates are ey |+-), |-+) and foe) which we 
number as ata ee 2y. | 3) and |4) respectively. The symbol | -+) re- 
presents the state where is of the nucleus X (whose spectrum we are 
determining) is -1/2 and in of the mamiike- nucleus A is +1/2. Using 


equation (40) we have 


Ca@ae ius oN Tone O Na /nr +R Gr 4 R Oo 

M2812) 7 oe eae 1212°12 * "1234°34 (44) 

= ha tug i Da J+ R +R 

Coder saa | Sana ny gee) haaga ade. 5412012 = 
Here Hb is (i/H 159, while Wi> = Wy J yy and Wes WytIayt™. Also 

= —baj 74 Se tee esi H iow 

O55 CE ve and abe UFEOred ais except Hoy or Has We used 

: a p : : { 
R B08 oR since oan if a=8. We have only considered the time deri- 

* e 6 

vatives of 19 and Oz, Since a is equal -to oer Also C14 and O59: 
which are respectively proportional to bs and HL can be ignored 


for low amplitudes of the rf field. Finally, we do not need to consi- 
der O14 and O53 Since the rf field does not cause transitions between 
states |1) and |4) and between |2) and |3). The differences between 
equations (44) and those of Smith's’ for this. case occur, in the. re- 
laxation terms, the R's. We recalculate these terms and then alter 
Smith's results for the lineshape function by taking into account. these 
differences. 


It follows from equation (41) that 


20 2 Z D has 
= -Yy-k We Ee 2yocy (al be Te ke, 45 
Ragas ~ ~3Y; {(8|Z; |s)+(a|1;|a)} + 215 (a jg CBIT; 1B)K. (45) 


: 3 
The first part of the right hand side of (45) is equal to - (3/2) k (y+ 


Y) which equals - (3/2) (Wy+W) where Ky is defined as W,, the spin- 


Ve 
MDs. : 
lattice relaxation rate for the X nucleus and RY is defined as ae 


the spin-lattice relaxation rate of the A nucleus. The second part 


om ow mst7 sana 


t 
¥ . 


anrelt 3) f vt 


he. 
bf SES 


of (45) is equal to (W-W)/2 with the result that 


Ree We owen (46) 


Meng 7 4 X 


By rearranging (41) for the case where a, a, 8 and 8 are all unequal 


we find 
D4 r] 4 . y] 4 
Re yy = Ly eki (ald layed, lB). (all aye ie. eye 
ene Me ee ee arco] aaleyce s 16 7) 
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C) Two Unlike Nuclei (and No Like Nuclei) Affecting the Resonance 


For this case there are three nuclei: the nucleus X (whose 
spectrum we are calculating), and the two unlike nuclei A and A. The 
prime indicates that the two A‘type nuclei are differently coupled to 
the X nucleus and the spectrum is called the XAA’ spectrum. In his 
calculations of the XAA’ spectrum, Smith (1972) assumed that Wy = Wa 
and that the A and A nuclei are not coupled. Since we wish to re- 
move These yresctrictyons 1t)is necessary fo recalculate the XAA’ spec- 
trum jab initio, — 

The eigenstates for this case (and case D in which we con- 
sider the AX resonance) are given in table 1(a). In table 1{b) a 
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Table 1 The Eigenstates and Transition Energies for the 


XAA’ and ARX Spectra 


Table 1(a) gives the eigenstates and Table 1(b) gives the transi- 
tion energies im units of h. The symbols D)vand 2a are dered 
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er EY is nearer Wy than Wy (otherwise the transition belongs to the 
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AAX spectrum). The frequencies Wy and w, are assumed to be far enough 
apart that the probability of the rf field flipping an A nucleus is neg- 


negligible near Wy. For the XAM spectrum we need only consider the 
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zante derivatives of o 36° "45? "46 and Or. The time deri- 
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vatives of the other values of ee are not considered because the 
/ 
transitions i<>j are part of the AAX spectrum or for reasons simi- 
/ 
lar to those given in part B. For the XAA spectrum we have 
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Since the eigenstates are orthogonal. Finally 
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Letting Gas = ae eS be the solution, we obtain a set 


of linear equations, which in matrix form is 
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where A, r and P are given in table 2(a). To calculate the spectrum 
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D) One Unlike Nucleus and One Like Nucleus Affecting the Resonance 


The eigenstates for this case are identical with the ei- 
genstates for the XAA spectrum and are given in table l(a). The 
transitions and the corresponding frequencies are given in table 
1(b). We shall merely quote the result in matrix form since the 
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method is similar to that used to find the XAA spectrum. We find 
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Table 2 Matrices and Vectors for the AAX and XAM Spectra 


The symbols S and C are defined as sin2a and cos2a respectively. 


/ 
Table 2(a) gives A, P and r for the XAA spectrum, Here 
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x(w)H, =-y fir-P. (66) 


where r is the solution of the matrix equation 
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Here A and P are given in table 2(b) and Z, = yeh’, /S2KT. 
It should be noted that (66) "gives the spectrum for both 
the A and A’nuclei, As we Require puie speccrulm tor ele eA mucous. 
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II(v) Theoretical Models of the Ruderman-Kittel Interaction in 


8-tin 


In this section we show how the theoretical Sn!!’ NMR line- 
shape in ®-tin is synthesized. We first consider the spectrum of 


each Sn!!7 


nucleus taking anor CCOUnt@,. InsOtaveas possible, the: inter- 
actions with other magnetic nuclei Yocated.in the first, second and 
third nearest-neighbour shells (since the Ruderman-Kittel interaction 
is believed to be of moderately short range, it is possible to neglect 
more distant magnetic nuclei or to treat their effect in a more ap- 
proximate fashion). We therefore need only consider the more common 


117 nucleus. Be- 


configurations or small groups in which we find a Sn 
Cause on Lheccomplexityv Orme newcal culation... We) lave SO tar aes able 
fOacerive expressions only: forescroups Of one; ‘two vand three nuclei 
(apart from the trivial case where there are four or more nuclei which 
ane all identical), “has 1ssnot a severewlamitation since the i1s0- 
topic abundance of the magnetic nuclei Snti°, Sn!!7 and Sn!!? is so 
smallp (see section: 11 Gy ))\thatethesprobabiiity.o: finding more than 
two magnetic nuclei in the inner three shells is small (according to 
Smith (1972), the probability is 19.7%). 

We consider three theoretical models for the Ruderman-Kittel 
interaction in §-tin. For each model we take different relative values 
of the strength of the first, second and third nearest neighbour in- 


teractions. However, before discussing each model individually, we 
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that both the x’ and A nuclei interact with the X nucleus but not with 
each otner. This contrasts with the notation XXA ingwhich there is in 
general an interaction between all three nuclei. We similarly dis- 
tinguish between the AXM and the XAA’ cases. 

The, first. model#that .wesshaidaconsirder 1s the nearest neigh- 
bour model, Here, we only consider interactions between nearest neigh- 
bour magnetic nuclei (each site has four nearest neighbours). We con- 
Sider the following basic configurations: X, XA, XAA, AXA, XXA and 
XA, For this model and the second nearest neighbour model (to be 
described) there is an interaction between the x" and A onuCleae Out 
not between the X and A nuclei for the Xx configuration. Similar- 
ly there is an interaction between the A and A nuclei but not between 
the X and A nuclei for the XAA’ configuration. Some of the other con- 
figurations have resonances which are equivalent to one of the basic 
configurations (for example, the Xx" configuration has a resonance 
which is identical to the X configuration, and the XAX configuration 
has the same spectrum as the XA configuration if there is no interac- 
tion between the X and X nuclei). The probabilities of nuclei having 
a resonance equivalent to one of the basic configurations' resonances 
are given in table 3 under the heading 'exact'. For this model, 
86.67% of the Sn!’ nuclei have a resonance equivalent to the reso- 
nancé of one of the basic configurations. For the remaining 13.33% 
of the Sn!!7 nuclei, we have approximated their resonance by the reso- 
nance of the basic configuration which it most resembles. The prob- 
abilities of these approximate resonances are given under the heading 


'approx.' and the sum of the exact and approximate probabilities are 
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given under the heading ‘total’ in table 3. 

The second nearest neighbour model is similar to the nearest 
neighbour model except that only second nearest neighbour interactions 
are considered. All other Ruderman-Kittel interactions are assumed to 
be negligible. As there are only two second nearest neighbour sites 
for each nuclear Site, the probabilities for the configurations list- 
ed in table 3 differ from those for the nearest neighbour model. 

The third model is the Ruderman-Kittel model. Here we use 
equation (34) to obtain the relative values of J. for each nucleus 
i and 7. Using the information en in the caption of figure 2, we 
obtain relative values of 1.0, 0.478, -0603 and 0.145 for Jas for 
first, second, third ‘and fourth nearest neighbour nuclei.‘ We see that 
the fourth nearest neighbour interaction is relatively small, thereby 
justifying our decision to treat explicitly only those magnetic nu- 
eller in the zoner chree “shells wiach consist O01 <0urmearest. Neien- 
bour, two second nearest neighbour and four third nearest neighbour sites. 

The probabilities of occurence of X, XA, XA’ and XXA con- 
figurations are given in table V of Smith's thesis and will not be 
repeated here (in arriving at these probabilities, Smith ignored the 
snits isotope because of its relatively small abundance). However, it 
should be noted that in 26.7% of the configurations involving three 
magnetic nuclei, the two magnetic neighbours are either first, second 
or third nearest neighbours of each other. We have taken such ‘cross! 
interactions into account though they were ignored by Smith. 

Finally, we take into account in an approximate fashion the ef- 


fect of magnetic neighbours outside the inner three shells. , The most 
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Figure 2 Structure of 8-tin 


The point group symmetry for f-tin is 4 2 m, with a = 5.831 A 
and c = 3,181 A at room temperature (Pearson, 1967). The num- 
bered nuclei are examples of the nearest neighbours to uekene 
A, the number n indicating the nth nearest neighbour. The num- 
ber of first, second, third and fourth nearese neighbour nuclei 
are four, two, four and four respectively. Using this struc- 
ture and the free electron model, we obtain a value of 2kp = 
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important effect of such nuclei is the broadening produced by the un- 
like nuclei. However their effect is approximately equivalent to an 
apparent increase in the relaxation rate of the X nucleus. This point 
is discussed in appendix IV. 

It will be seen that we have ignored interactions between the 
more distant magnetic nuclei and those located in the inner three 
shells. It can be shown that this is quite a good approximation (D. 
G. Hughes, private communication). The effect of mutual interactions 
DetWeen mucici "Outsiae theriinersenpeessnel ls are, expected to pe ¢- 
quivalent to a small reduction in the average spin-lattice relaxation 
"fate of “such nuclei. “Such"€ffects are taken into account by the in- 
troduceion ot 4a a1tting parameter co as.aescraped in appendix, LV. 

The theoretical spectra which are required to synthesize the 
Sn1l7 lineshape are found using the results of section II{iv). To 
find the X and XA lineshapes, we use equations (43) and (49) res- 
Pectively, the XAA’ lineshape is given by equation (65) in conjunc- 
tion with equation (63). The XXA lineshape is given by equation (66) 
together with (67). The solution of equation(65) is given in appendix 


II while the solution of (66) is given in appendix III, 
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ITI APPARATUS AND EXPERIMENTAL PROCEDURE 


In order to study the Ruderman-Kittel interaction, samples 
of tin were rapidly rotated about an axis inclined at the magic angle 
to Hi While this could be done by spinning the sample about a ver- 
tical axis and tilting the magnet until the magic angle condition is 
Satisfied, we chose to spin the sample about an axis inclined to the 
Wereical’ soO,that the macnetscouldibeskept Heunanenl, 

The high speed rotation was achieved using an ae turbine 
Similar to that pioneered by Beams (1937).-. In this system air under 
pressure passes. through fine holes or jets drilled at. a suitable angle 
in a conical stator and impinges upon flutings machined into a coni- 
ealerocor, wlhese jiets.ot air Unt tythesrotor orn the Stator and (start 
it—spinnine.; The rotor does not fly outsof the stator, but: rides u- 
pon a thin cushion of air just above the stator surface, as can be 
understood by the application of Bernoulli's principle. According 
fon UiaseOrinca pIlcothe pressures on tne minders idesOr) Che TOLOn 1s sre 
lated to atmospheric pressure by the relation 


+ aii =. P (1) 


faowex atmos 


where p is the density of air and v is the VelCcity, oF tie ais .iu 


contact with the underside of the rotor. Since Citeuee is smaller 


than P there is a downward force on the rotor which tends to 


atmos’ 
keep the rotor in the stator. Indeed, this force even enables the 


rotor to spin upside down. 


Since the turbine had to be operated within the rf coils of 
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the Varian NMR probe, the stator and rotor were made of nylon. Most 
of the rotors were made out of Delrin (Ertacetyl H) which has better 
mechanical properties than ordinary nylon (much larger moduli of' e- 
lasticity and about 20% greater tensile strength). The stator was 
constructed so that the axis of rotation of the rotor made an angle of 
40 degrees to the cylindrical stator support which just fitted within 
the Varian probe. ‘Because the :stator Support was perpendicular to Hy» 
the axis could make any angle between 50 and 90 degrees with H,, mere- 
ly by rotating the whole stator assembly. The magic angle (54°447) 
could therefore be obtained by suitably orientating the stator. 

Piesstator previously Vescripecd DY (omni (louc)) Was diseq. 1x, 
the present measurements. However, our rotors differed substantially 
from those used by Smith, and the design is shown in figure 3. They 
were hollowed out so as to achieve a maximum sample volume consistent 
with stable rotation within the Varian probe (inner diameter 1.7 cm). 
ihe sample volume of the rotors ranged between, 0.23 and),0.25 a and 
was 2 or 3 times as large as those used by Smith. Another point should 
should be made about our design. The center of the plane through x'-x 
and perpendicular to the axis of symmetry of the rotor; a, 1s approxi- 
mately in the same position during rotation. The center of gravity, 
@. Opviousiy lies alone athe axis, of rotation se AS a cesult,= che axis 
of rotation goes through points g and a to a good approximation. Thus, 
large perturbations result if the angle between the line g-a and the 
line a-b is large. This can occur for slight asymmetry of the rotor 
if g is very close to a. In our design, g was located as far above 


the x-x' plane as possible. We attribute the improved stability of 


Ass 


i. 
os Ge - i toh 


° 
We oly : for 69 
i 
" ) 1 ifie 
ins. yv a 7 4 
yi ' ‘ 
L 
¥,,;9 ; 
. 4 é % 
Ye iy (A #D 
; al j 
{ il) ” 3 
id 
| ; 4 , 
Oe ai . 4A 
‘ 2 ‘ 
( ; + ‘ . 
{ 
a ‘ y i 
- > 
fe : | rai% ej » i? } 
r 
, 7 6 i A 
i . A = 
ids int , t-05 LP 4 
: - ; + 
‘x Ago ree te j 
> 
a ; an = ol 
=F, : € 2 2 08OT i a 


(itlvesz & ZA 


sari ‘ re bye 


U 0 


eed hooy 


ree Bare I } 
: in sont ads ho Yet ome 
pe 


i 
ae 


| me a 

‘tel tw irene > wivisd Ye’ 

' _ i a i 
ty 


spre stant) poetvit 


Y tS ye 


«MOLTEIG 
_ 


: ‘em bee s-y anil edz woowred af 1 oddt 2 Slurees 


ent xt 


all phsed nee 


egret ine “th 


¥ 


pee neat a sci ten > anaion 


rsa wiieney + HOE tins ds dada neois 29 a 


5 i 
Sisnts 


LP | 
> { “i: eF 
A oe 2 Ue 


_ rd 
© : AD = ee iz Fi + hod 
1 
. , 4 . 
* h, « (se a * > a 
a. : a ds) a? Oo a We Cuore 
S* 
7 yy = \ ats 
: j » A« 
~ ofam bludp 
»» a eet i” ‘ 
( | & ant i ‘ 
—Ay 7 
. - mq 
» 4 J : = 
2 f r 714. + Lq 3¢ ¢ > 
Gg 
f ; per? cy t 
rt — ‘ 
‘ ’ e ri 
i ? cng my 
-— 
WIS Oe i tT } ¥ o4¢ 
— 4 ee, 
- Wis 25 Det ty = Yo 
; r 4 ) 
ved ‘ —~ - > use 
ss TiIsguve TH ry ry a3 5 2 LORS Teq 
' i : 
: . . y r’ * i b 
mitatot anim agitican omaz o1i3 i 


to sins ona snaltc eoll eteuotvd ” 


- 
oF £ oo sinter guerre 


rsag noisere 
= 


io 


Pit 7 al a an 


104, ory e ‘omuget hy ee 
| A ve “nm fe ieeh ; is mi i ; 7 
S 007s RY i ee dati 


‘Mion bk i ae of a eral oe ee er 
Cie i ; ays M oe i oh as) nis ; Rien | ay 


Aas ot a how si h 8 A al y 4 Lah! yy, 


ai nk. orf aida Atty: 8 


1 rat a a a nv) : , i Loos Bar Hai | 
ee 9 saiaDee wy 20! nok ar tated 7 ne wie nnitey eu ae a) tok onan 


a. ‘on 4 are ; me: pay Y a iF aed 1 


(i! ‘ 
- iy 


Ser 7. fot q nore) e 2 : a J to ah ag ia (3, ‘ vN +. axes t ha par | \ Pik ih i 


ea 


ie bee Y ok Mt A | i bi cy ) mae fing i iy Ne et 


Lal ; 
i : 
Py ie | 
“a 
Piginh i } 
iRP i Mg 1 hi Me 
Vy f i 

4 j 
i 

oy fa i 

af 
an! : 1) ul 


Pek ee) 


———— ERO ye eng S 
—s Zi . 4 7 _ a i 
o = s = 


ib ie ay i av | ay 
" ‘\ a | ot i? ‘ he ; 
ne ey Panay 
ioe | : . i ioe 
ci oe oa ah 
ni, - a aa tae ay 
4 a ae ot ie ‘ Te var 0 


ee | : " . a i 


a y 
a b, 


oy 
iy an 


Figure 3 § The Rotor ‘and Stator 


Here a, b, x, x and g are defined in the main body of the thesis; 
Jets of air pass througha from id; ¢ 1s the stabilizing {relief} 

jet (Smith, 11972); f£ gives the position of the walls of the anter- 
ior of the Varian probe; g’ gives the, position ef the receiver coils; 
heand a represent “the two parts of the stator/wnhich are glued to-— 


gether (omith,)£9/2)< 
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our rotors largely to this modification. We also found that a more 
stable rotation of the rotors was achieved when the maximum diameter 
of the rotor was less than the diameter of the top of the conical sec- 
tion of the stator. Both eight and twelve flute rotors were used. 
However, little consistent difference between the two types was ob- 
served. 

In order to achieve full penetration of the rf field into the 
metal, the samples consisted of powder set in epoxy glue. The tin was 
purchased from) Fisher,ScientifieyCompany. «Fhe purity is beli.evedpto 
bein the range 99.9 to 99.99%... However attempts to confirm this 
with Fisher Scientific have so far failed. In order to accurately 
determine the magic angle setting, a small amount of aluminum powder 
(roughly 10% by volume) was added to the tin powder and the magic 
angle was obtained using the Al*’ resonance. 

The advantages in using the Al*’ resonance to find the magic 
axis are: 


(a) The spin-lattice relaxation time, T,, for Al?’ is approximate- 


1? 
ly 6.3 milliseconds at room temperature so that the linewidth 
(defined as the interval between points of maximum and minimum 
slope of the absorption line) associated with the Ty broadening is 
30 Hz (compared with more than 1 kHz for Sn!17 and Sn!19). 

(b). The natural abundance of A127 is 100% which means that the 
Ruderman-Kittel interaction should not broaden the Al*’ resonance 
in pure aluminum (at least intprinciple)). 

(c) The peak to peak derivative linewidth in a stationary 


sample of aluminum is 9.4 kHz (Gutowsky and McGarvey, 1952), where- 
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as the corresponding linewidth for Sn and Sn Imspune timiis 

3 kHz or more (the precise value depends on the field strength be- 
cause of the anisotropic Knight Shift). 

(d) For equal numbers of nuclei, the NMR signal strength for AIS! 

at constant frequency is more than 8 times as strong as that of 


117 whereas 


Snl!17_ Moreover only 7.67% of natural tin nuclei areySim 
100% of aluminum nuclei are Al@’ (NMR Tables, Varian Associates, 
5th edition, 1965). 

It should be noted that there is a possibility of a nuclear 
quadrupole interaction in aluminum on account of the spin number 5/2 
of Al?7 . Lattice defects or strains mayveivesrise tornon-zerolelec- 
Ericarleldsgradiéntseatethernuclearvsitessw) However; sthrsyquadrupole 
broadening should average out provided the sample is spun sufficient- 
ly rapidly. 

In making the samples, sufficient epoxy was used to avoid air 
bubblessin' the mixture, since these would seriously unbalance the. ro- 
OTs? 

The equipment used to supply the air to the turbine has been 
described by Smith (1972) and will not be described here. Standard 
techniques were used to measure the speed of the rotors (Henriot and 
Huguenard, 1925). 

The NMR signals were observed in the absorption mode using 
a commercial Varian wide-line spectrometer, model VF16B. It was found 
that the noise output of ae spectrometer increased when the rotor 
was spinning. The excess noise was found to extend over a wide fre- 


quency range so that it could not easily be filtered out. Its magni- 
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tude was found to be proportional to the rf voltage across the trans- 
miGter, coil, indicating that the source of; the mrben was a microphonic 
vibration of the NMR coil system which modulated the coupling between 
the transmitter and receiver coils. It was found that placing cotton 
wool near the top of the Varian probe (where the sample was introduced 
into the probe) reduced the microphonic effect, presumably by damping 
out the acoustic vibrations. Also, masking tape was placed along the 
inner walls, of the’ Varian probe: (except over the receiver corl) to 
further reduce incrophonics. “AS asprecaution, we: reduced. the audio 
bandwidth of the spectrometer immediately following the detector from 
about 10 kHz to about 500 Hz to prevent overloading due to the large 
Signal present at the frequency of rotation (approximately 5 kHz). 

The NMR signals were recorded as the first derivative of the 
absorption by using a sinusoidal field modulation at 37 Hz in conjunc- 
tion with lock-in detection. The output of the lock-in detector was 
fed to a 1024-channel digital signal averager, Fabritek model 1062. 

A block diagram of the equipment is shown in figure 4. 

Synchronization of the magnetic field sweep and the internal 
sweep of the Fabritek signal averager was achieved by triggering the 
signal ayerager with a signal from a photodiode which was turned on by 
light reflected from a-piece of aluminum foil attached to the rotating 
field-sweep potentiometer of the Varian Fieldial regulator. Repeated 
sweeps through the resonance made during a 12 hour period were stored 
in one half of the memory of the signal averager. This procedure en- 


abled us to cancel any slight asymmetry in the signal caused by the 
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lock-in detector, by switching the polarity of the signal fed into the 
signal averager from the lock-in detector after each twelve hour run. 

The results of each 12 hour run were then transferred (with appropriate 
polarity) into the other half of the signal averager's memory. After the 
the entire series of runs was completed, the resonance curve was read 

out onto an x-y recorder and a theoretical Mneshape was fitted to it 


in the manner described in the next section. 
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IV RESULTS AND DISCUSSION 


IV(i) Experimental Sn!17 Resonance in 8-tin 


The derivative of the Sn!!7 NMR absorption signal in a station- 
ary sample of §-tin at room temperature (22-C} 4g shown in figure 5. 
In obtaining this resonance, H, was 7000 gauss and the magnitude of 
the magnetic field modulation was 2.32 gauss. The asymmetry of the re- 
sOoliance 1S Characteristic of thefanisotropace Knight Shift inja poly- 
crystalline sample in which the nuclei are located at axially symmet- 
ric sites. The rf transmitter voltage was adjusted so that satura- 
tion effects (Bloembergen, Purcell and Pound, 1948) were negligible. 

The effect of rotating the sample at the magic angle with an 
angular speed of 5.4 kHz is shown in figure 6. The resonance is the 
result of summing 25 twelve-nouriruns. The resonance 1s quite sym-. 
metric showing that the anisotropic Knight Shift has been averaged 
out. Rotational sidebands, separated from the main resonance by 
5.4 kHz, are clearly visible. In order to obtain a reasonably faith- 
ful reproduction of the absorption derivative, the modulation ampli- 
tude was. reduced to 0.32 gauss in obtaining this resonance. Also, 
the rf transmitter voltage was reduced to such a level that the loss 
of peak-to-peak intensity due to saturation amounted to less than 3%. 
Both these adjustments reduced the signal intensity. 

Since we wish to compare the lineshape of the experimental 


resonance with theory, it is necessary to correct for the residual 
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instrumental distortion caused by the finite sparen amplitude and 
the rf transmitter voltage. The corrections which are applied to the 
theoretical rather than the experimental lineshape have been described 
by Smith (1972). The value of the "saturation" parameter p, defined 
in Appendix IV of Smith's thesis, was .025 for the resonance shown in 
figure, 6. 

Finally, the experimental resonance was corrected by subtract- 
ing the rotational sidebands. These were assumed to be of the same 
shape as the main resonance, and their amplitude was estimated by vis- 
ual examination of the resonance in figure 6. The experimental reso- 
nance, corrected by removal of the sidebands, is also shown in this 


figure. This is the resonance which will be compared with theory . 


IV(ii) Procedure For Fitting the Theoretical Lineshape to the 


Experimental Resonance 


dhe Veross-over sora centrals frequency , andthe; baselianeiof 
the: experimental, resonance were first determined. The-'cross-over" 
point can be found quite accurately on account of the steep slope of 
the absorption derivative near the center, and it was assumed that 
this,point coincadesswith phesicenter ef ithe: theoretical dineshape. 
Next, the experimental resonance was digitized by measuring the in- 
fensa ty f. (relative to the previously chosen baseline) at a series 
of equally spaced frequencies v; on either side of the cross-over 


point. A total of 61 points were taken consisting of the origin 


(cross-over) and 30 points on each side. Since the interval between 
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points yao aaa was 141 Hz, the fitting region extended over 8.46 kHz. 
To take account of the different intensity of the theoretical 

and experimental lineshapes, and to allow for possible error in the 

baseline determination, the theoretical lineshape intensity at fre- 


quencies v; was expressed in terms of the normalized theoretical line- 


shape function g(v) in the form 


h(v,) = a, + a,8(v,). (1) 


The quantities ay and a, are theretore, ire (baseline correctiom and 
normalizing factor respectively. 
The fitting was done using the least Squares criterion that 


a ol bea aieny (2) 


is a minimum with respect to the fitting parameters. The quantity 9; 
is the standard deviation of the ith data point as estimated from the 
Signal-to-noise ratio of the experimental resonance. 

The normalized theoretical lineshape function g(v) is a func- 
Fouad uCOtal some N|2 5 


tion of the fitting parameters Ons O a 


4? N’ 
while the lineshape function h(v) is a function of N parameters. For 
example, in the nearest and second nearest neighbour models, N is 3 and 
O is the. value of the J coupling between nearest and second nearest 
neighbours respectively. In the Ruderman-Kittel model, N is 4 since 
there is an extra fitting parameter Oy which takes account of interac- 
GIONS sWi Gh dis cant-nuclei.. In ‘this case, oi is equivalent to the para- 
meter € defined in section II(iv). 


For the fitting procedure, we require the values of BEY 5G ne 


However, since the calculations of g(v;) in general involve the inver- 
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Sion of an 8x8 matrix (the matrix A in table 2(b)), this has to be done 
numerically for each value of i and for each value of a , and it is 
J 


not possible to express g(v) as a simple function of the a.'s. We 


’ therefore approximate gv.) in the form of the Taylor expansion 


N N 
Oy) corks ta ep oOe + ye C. «16a. da, (3) 
1 af 423 oh ee jk=3 1J J 

where 

Qe sae) ON. 4 

J jo J 4) 

and oO 1S a-suitably Chosen value, or oe The (CoOerErcrents A;; ie 
J 


and oe were found for each value of i for the nearest neighbour mod- 
1) 

el (where N = 3) by fitting nine different explicitly calculated 

theoretical curves to equation (3). The nine theoretical curves cor- 


Poe OL. Fo OOM rn re aeRO Krzr Os 


responded to val £ =p! 
p e o values o a, | jo 


Deingucdual tO 9.0 “Ki, DINcemune COCLLICICNLS are overdetermined, 
cern Bese’ Value Was foumidsusing a Lease squanes tit. “AS ‘a "check 

On the ValtdTEveor  (a)) Over vsuchwal extended, ranve,of Os the nine © 
theoretical curves were reconstructed using the expansion. They were 
found to be in excellent agreement with the original explicitly cal- 
culated curves. The above procedure was then repeated for the second 
nearest neighbour model, For the Ruderman-Kittel model where N=4, a 
total of 18 theoretical curves were used to calculate the coefficients 
Ais ae and C.. 


i ye 


IV(iii)- Results 


The fits of the experimental resonance of Sn!17 for the nearest 
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neighbour model, the second nearest neighbour model and the Ruderman- 
Kittel model are shown in figures 7, 8 and 9 respectively. For these 
fits we used a value of T, = 169.5 microseconds. This was obtained using 
the measured value of 155 + 3 microseconds for Sn!12 in g-tin at room 
temperature: (295°K) (Dickson, 1969), and using the fact that Ty) is in- 


versely proportional to the square of the gyromagnetic ratio and that 


Nema an 17 = 1.046535 + '.000003 (Smith, 1972). The value of oO. 
is taken to be .04 intensity units of figure 6 except for the cross- 


over frequency and the three frequency values on either side of it. 


POretnese seven frequency values, o- 


j 1s taken to be .32 intensity un- 


its . The fits for the nearest and second nearest neighbour models 
are poor with e being equal to 5860 and 12,750 respectively. How- 
ever, for the Ruderman-Kittel model the fit was good, the value of 
ce being equal to 127. For this model the value of J, where J is 


the value of Jee i and & being a pair of nearest neighbouring Snll?7 


andisn++? 


nuclei, 1S 35: U4 ys 0omkiz. (thissiimplies sthat.J 71stequat °to 
2.91 + .06 kHz for two Sn!!" nuclei and 3.18 + .06 kHz for two Sn!19 
nuclei). The error is determined from the fit (Guest, 1961). The 
value-of € is .49 + .05, the erxor again determined from the f1t. 


We have not taken into account other sources of error; this is done 


in section IV(iv). 
IV(iv) Discussion 


It is clear from figures 7 and 8 that the first and second 


nearest neighbour models are in very poor agreement with experiment. 
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It follows that neither the first nearest neighbour nor the second 
nearest neighbour interaction is dominant in 8-tin. It has been sug- 
gested (Sharma et al,1969) that the second nearest neighbour inter- 
action could well be dominant, on account of a hole surface in the 
fourth zone of the Fermi surface with fairly flat faces perpendi- 
cular «to the (o01] direction. According. to the results‘ot Roti et 
al (1966) this should give a strong interaction in the foo1] direc- 
tion, that is, in the direction of the second nearest neighbours. 
Qur=results show quite definitely that this is not ‘the:case. 

berCana De: Seen ei elCUre oi natetne acreenenusne tween sexpert — 
ment and the Ruderman-Kittel theoretical model is very good. The 
"goodness of fit' parameter x? was found to be 127 as opposed to the 
"expected' value of 56. However, such a discrepancy is by no means 
unreasonable since we estimated the magnitude of the noise by a vis- 
ual examination Of figure 6, am unreliable procedure.’ On the basis 
Or tne-results presented in tis, tuestseatnere.as nO reason, co .doubt 
the validity of the Ruderman-Kittel model as applied to §8-tin. 


it will be moted thatvour value OfsJ depends upon the values 


Tey Ne ES, 


of qT, which we have assumed for Sn and Sn in 6=tan., in order 
to check this point, we repeated the fitting procedure by assuming a 
2% error in the qT, value as given by Dickson (1969). However, the 
Value of J was found to be changed by only 0.05 kHz, so it is clear 


that uncertainty in T, is not of major importance in our determination 


At 
Of-os 


Another possible source of systematic error in our J value a- 


rises from the difficulty of estimating the amplitude of the rotation- 
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al sidebands. Any error in the sideband amplitude will obviously af- 
fect the corrected lineshape within the fitting region which extends 
4,2 kHz on either side of the cross-over point. We estimate that the 
error in J arising from this uncertainty is about 0.15 kHz. 

Another source. of systematic error in our value of J arises 
from the various approximations made in synthesizing the seheane tical. 
lineshape. We estimate this to be about +.25 kHz. By combining the 
various errors we estimate our value of J for the interaction between 
Sn!17 and Sn!!9 nuclei which are nearest neighbours to be 

Jue O04 ee Oe Sakiizes 

We note that the Ruderman-Kittel fit was really a two-para- 
mecenet Lt Since Uiexpatanetcnm earcorescne ince newetrnect wor distant 
nuclei was allowed to differ, from the value of unity. The fact. that 
& was found to be 0.49 +.05 indicates that the broadening caused by 
unlike magnetic nuclei outside the inner three shells is significant. 
We attribute the fact that € is less than unity to a combination of 
the following reasons: 

(a) in the derivation of equation s(5) gin appendix IV) 3it)was.as- 
sumed that the linewidth contributions due to each distant magne- 
tic nucleus add linearly... However, thissis only truesit the line- 
shapes associated with these interactions are Lorentzian (Hughes 
and MacDonald, 1961) and if they are independent. In actual fact, 
the lineshapes converge slightly faster than a Lorentzian function 
(D.G. Hughes, private communication), and one might therefore ex- 
pect the resulting linewidth to be somewhat smaller than the sum 


of the. individual Jinewidths, that.is. one would expect € to be 
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less than unity. 

b) Mutual interactions among the distant magnetic nuclei outside 
the inner three shells have been neglected. However, mutual spin 
flips associated with the Ji,.L, interaction will, if the nuclei are 
of like species, tend to reduce the broadening cause by these nu- 
clei, by a ‘motional narrowing' effect. This again would tend to 
LeEGucere. 

(c) No account has been taken of the experimental error in the T; 
value given by Dickson (1969). However, it is intuitively obvious 
that a 2% error in Dickson's value would have a large effect on &. 

Apart from the value of 1.89 + 0.09 kHz found by Smith (1972) 
for the strength of the Ruderman-Kittel interaction between nearest 
neagnbours in-p-tin, che following values appear ineethe Literature: 

22025057 kHz™(MeLachian, 1968) 
2.5 kHz (Karimov and Schegolev, 1961) 
Ani et OS SR KHz CAT loulmand= pel tours. 1969); 

Also Sharma et al (1969) concluded that the ‘exchange narrowing' they 

observed with a crystal of isotopically pure Sn!!2 was consistent 

with Alloul and Deltour's value. 

A major drawback of the various methods used by McLachlan, 

Karimov and Schegolev,:Alloul and Deltour, and Sharma et al, is their 


119 LeSO= 


reliance on the measured value of the second moment of the Sn 
nance in natural tin at low temperatures. The difficulty in measur- 
ing this second moment is illustrated by the fact that Karimov and 


2 
Schegolev found that the second moment is 1.2 + 0.3 (kHz) whereas 


2 
Alloul and Deltour obtained:a value of 2:5 +(@3 (kHz). 
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In addition, we feel that McLachlan's value of J is unrelia- 
ble since he assumed that the crystal structure of §-tin is hexagonal 
close packed and that only the nearest neighbour Ruderman-Kittel in- 
teraction was significant. 

Karimov and Schegolev, on the other hand, ignored the pseudo- 
dipolar interaction entirely. This appears to be an important omis- 
sion since McLachlan found that the pseudo-dipolar interaction was of 
the same order of magnitude as the dipolar interaction. 

Alloul and Deltour's data, obtained using a spin-echo method, 
is unfortunately by no means easy to interpret. We therefore feel 
that any beara between their value and ours may well be due to a 
combination of the approximations made in their theory, and any error 


in the measured value of the second moment of the Sn!19 


resonance. 
Finally, we ask why the value of J obtained by Smith should 

be so much smaller than our value. We believe that the main cause. 

was the poor signal-to-noise ratio of the resonances obtained by Smith, 

together with the restricted field sweep used by him (10 gauss com- 

pared with 25 gauss in our case). The rotational sidebands are bare- 

ly visible in Smith's resonances because of the poor signal-to-noise 

ratio, and they were not corrected for. The ‘wings' of the resonance 

therefore converge more rapidly in the fitting region because of the 

sidebands, thereby giving the small value of J. Moreover, Smith 

did not do a two parameter fit as we did. Rather, he had to assume 

a Value for the ee eeeree (proportional to Bi an'his notation)... “the 


value he assumed turned out to be larger than the one we found. This 


again would tend to give too small a value of J. Also, in synthe- 
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Sizing the theoretical lineshape Smith made several approximations 
which we have avoided, though it is difficult to tell whether these 
would be likely to give too high or too low a value of J. 

In order to confirm our value of the strength of the near- 
est neighbour Ruderman-Kittel interaction, we suggest examining the 
Sn!!19 resonance in the same way as was done for the Sn!!7 resonance 
in this thesis. Finally, it should be possible to computationally 
broaden the Sn!!’ resonance we obtained with the spinning sample, 
by the anisotropic Knight Shift and direct dipolar broadening, both 


A 


of which are known. A comparison with the Sn resonance of the 


stationary sample would then reveal whether the pseudo-dipolar inter- 


action produces significant broadening in §-tin. 


65 


ae | mh et 
S408? YuLrorye 
a ne 
’ ¥ 


i ’ A a S or 7 
oly wk oomenbsom €4 ¥ 2 
er ake _ » iy 
ni} -.esnedt eld ak 
7 i 7 nd a 


REFERENCES 


Abragam A., 1961, The Principles of Nuclear Magnetism, 
(Oxford: The University Press). 
APVouls Hey and Del tour’R., "1969, Phys. Rev. 183, 414. 
pHaecrsOnger. Wes, 1954.0). Phys. Sce- Japan 9,..516. 
MICRE Wr Rvs) DradoUry Ape atid. hades h~G.,. 1950, Nature 152 1659, 
Did rew PE shnernaera ried Lelie, eLoOoMmiO bec, PHYS elo je) los. 
NidreW SEER and Jenks) GJ. 1962, Proc.) Phys, 0c.) 60, 665- 
pnarew ERs ands Newang R. An foo, Proce Phys: 90c. 72,7959. 
Beaiismon Weslo) sede ADDIS Jriiy st. ctu 795". 
Brochets, HanseneWeW: and "Packard MeEt5 1946, Phys. Rev. 70,474. 
Bloembergen Nv 5oPurcell £:°M.%and*Pound’R.V~*,~ 19485 "Phys. “Rev. 
75, O79. 
Bloembergen N. and Rowland T.J-; 1953; Acta Met. 1, 731. 


Broembergen:N.”and* Rowland T?J=, 1955, Phys: Rev. 97, 1679: 


Davydov A.S., 1966, Quantum Mechanics, (Ann Arbor, Michigan: 

NEO Press). 
Decksonep M1000) Phys kev 31840224. 
Guest P., 1961, Numerical Methods of Curve Fitting, 

(Cambridge: University Press). 

Gutowsky H.S. and McGarvey, BR.) 1952) 2 Chem. Phys. 20,0 1472—° 
Henriot H. and Huguenard End 1925, Comptes Rendus 180, 1389. 
Hughes D.G. and MacDonald D.K.C., 1961, Proc. Phys. Soc. 78, 75. 
Karimov Y.S. and Schegolev I.F., 1961, JETP (English Translation) 


13, 908. 


66 


oe ho) abies heer am 
' (zemnt eters ie He iano , 

bee eat ok att | oer le theston b ba oy tuo tA 
ate 8 waiqal. 08 ‘euie . jeaer + 4.3 Hoes 
W2at .t8) warrar ‘ee * S, 4‘ weball bee vA a “|, a he 


90.,.08 <cek cevidld cert .teer | ub tanel bts Ws rt iach . ; 
oe ‘ . ' i ar) 
east. ; 7 > “7 
‘SER oY cee att) 5024 8287 .. A gobs bee 1d een 


tet A ew J CqGk Teed WL eae 


< ey = : +1 = 
; ‘ : ne - : 
PSS. OF VOT edd .eehl 0M bintoh? San. 8 aseael yt dsolt =e 


WOR 2 ,SEOY , A.4 feo bas Med lkesan ‘av tag7 sUaubo [E 7 : 
“7 ° : Bp la 
; ; — ove AY Ve SB . : 7 
TEN. ft .Del aaah ,c2e! ,.L.7 baalwah See ff aug isdwed tS A a 
ed o > , _ 
: $ J 
eta SS VOR ,ayet 9200 ,...7 hootved bre of asevedmsole ‘ 
! 4 . ‘ - 7 : E 
taegitsie ery intA).29inclooM mincaul .daer , BA vobyynd if, 
y 7 - 
| , fener? om, Soa ; 
Alia Al | Ly OR ME 9 RAY RHEL", Hae aio | ; 
| ’ 7 an ny | ; : if : 7 - 
; cite Sn ek ee Ui gnosis oven ba atog sy Leojrreigait fees =n Balk ; 
ie 6 7 apa ; 


hn ee oe 1 eleeer saterovien sanhsdimeD), \ a 
Rane - a vba) aan a , 
Lore, ot; see Le eae ea: 


— 
vie ey fi 
a : ' 
- _ sya ne ae 
ae Lida he erect 


67 


Kessemeier H. and Norberg E.R., 1967, Phys. Rev. 155, 321. 

Knight W.D., 1949, Phys. Rev. hee As 

Mahanti S.D. and Das T.P,, 1968, Phys. Rev. 170, 426; 

McLachlan L,A,, 1968, Can. J. Phys. 46, 871. 

Pave Go. 547 1250, oOddid olate Physics 2 yeuNew Yorks AcademiesPress) , 

| pre 

Pearson W.B., 1967, A Handbook of Lattice Spacings and Structures 
of Metals and Alioys*2,)\(lerento:) Pergamon=Press)} pe 25: 

Popie J A; Schneider W.G. and Bernstein Hod fF lO59 miter Resolutrom 
Nuclear Magnetic Resonance, (Toronto: McGraw-Hill). 

Purcell E.M., Torrey H.C. and Pound R.V., 1946, Phys. Rev. 69, 37. 

Kedivelarn G75 1957; PiBM Journal ely e19. 

RecuMine ecoleer Heo. eanG) KADPAN Sens, “1900s aL NVSuanev mito e olon 

Ruderman “MVAS "and "Ketter ' Cy, 91954, "Phys. +Rever Jo, 99. 

scuwind “A. 52>) 190/>5eAnn. Pnys@k wy, 922% 

Sharma S.N., Williams D.L. and Schone H.E., 1969, Phys Rev. 188, 662. 

Stichter*€7P. 5719635, ‘Principles ot Magnetic Resonance; * (New* York: 
Harper and Row). 

Satta wih. Lose ee hepee ests, Unvversicy, 01 Al berear 

Van Vieck J5H., 194877 Phys. kev. 74). 1t03er 


Winter J., 1971, Magnetic Resonance in Metals, (Oxford: Clarendon 


Press), pp: O/. 


oer 
~) ' “ Ean ye : ot) e . 7 ; ae @3 ; bates 


‘ ! yshi¢ 
cea 
a5 ‘ 
7 Ly rst } wor 
‘ ; J P24 t 4 ies 
ee ee oy 
i 
. ; : a list 
z 
. r ™‘~ 5 
‘ ‘ é . 
Pe 
a CsA 
$ «08! 
/ J : 4 7 ) i « we tj : . ; -— © : 
: 
ol s i 
io wai") ‘ i 2 2atvae 2 2 Saad | ~Lae 
Fie: 
> bt 
Sieg 
t * «2 
Bt pedis o VYile€tavs U ,28eadt Pe P 
® 
. i 
oaTT a Ww 2 any it ar -" ah | 
~ereld at .WoaR ea at ® ot cat. 1. doe — 
. ; - [ 
‘ Ue 
J 
ae - aM ; sid : = Sr. o 
| mobaowsdy sbeotino) etarel Al sdetnsest sedans pbedtet Tee 
® ; 7 


7 7 5! i aay Pea Sd cee 


pe 


va res 
iat ue ‘= : vl . Ps 
5 i. 7 


i - 
©) xs 


68 


: 2 
Appendix I The Time-averaged Values of Co 


. ; : : 2 
In this appendix we find the time-averaged values of Co where 


CoH is the direction cosine for the nth principal axis of a symmetric, 
second rank tensor T, with respect to the z axis (the steady magne- 
mic eee dudirection)) iicould reptesent the Knight onirt tereee of 
8-tin for example... Suppose ‘that the principal axes of T are x yand 


M4 
z. AS we want to find the time variation of Oh (in order to obtain 
TY 5 


] YW W We 44l 
the time averaged value of oe ), we let x, y and z rotate about an a- 
n 


' 
xis z with angular speed (2 “Thewane le between Zz and’ zis Se By de- 


termining the vectors x. aoa oe terms of the fixed laboratory 
coordinate system xyz, we can evaluate the time variation of Co toince 
the above vectors are unit vectors and thus Oe is just the z compo- 
nent Or Che Nth princi palwaxis ans the xy zycoordinate system). incor 
der to do this we introduce another coordinate system xyz where 2 has 
been previously introduced and the yz plane Corresponds utor chesyz, plane. 
Death principal axis sis thes vector (sinx_sin(at+y_),sinx cos(at+v_), 


s¢t : : 
cosx.) in terms of the xyz coordinate system. Here x is the angle 
n n 


between the nth principal axis and z, and » is the angle which is 
n 
needed to fix the orientation of the nth principal axis at t = 0. A 
4 o¢7 : A 
wector Ain) the xyz coordinate system 1S the vector A in the xyz coor- 


/ 
dinate system where A = R,A and 
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Thus the nth PrincipdalFaxts 1s tne vector 


(sinx sin(2t+) ), cosa Sinx cos (t+) )-sina COSX, 


& 5 


Sina _siny cos(Nt+ +COSa_COS 
asinx, ( ‘ie O. x? 


in terms of the xyz coordinate system. Therefore 


Cr = sine Siny COS (etry ) “+ COSa COSY 
on 3 n n 3 n 


and 


E Si é tne /2e+ wee aie 
= S17) i a 
3n aS Mp 3 ue : 

2 
When a. satisfied the magic angle condition, sin a =e 2/15" iTCeCOS Be = 


Z 
L/S withethe result that C 
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Appendix II: The Absorption Lineshape for the XAA Resonance 


In this section we include the expression for the absorption 
lineshape X(w) of the XAN resonance. This was obtained by solving (63) 
for r and then substituting the solution into (65). The x {w) was 
determined from x(w) by noting that fw) is equal to the imagin- 
aryepartvort 71) taetTher soluta ones 

x'(w) = (ovat /L6KTW, ) [B{ATL+g (F-2D) -2A-4G} + co” {B(2A+4G+g (2D-F)) + 


4 4 
J (EF-2ED-4H-2AK)+2Ab } - 2A(be) | 


* 2 4 4 
B {JC+4A-4gF+8G}+ c {B(4gF-8G-4A)+J (4AK+8H-4EF)-4Ab }+ 4A(bc) 


where 
2 Digs es hvak 
A= (Ag-e) + 4hd Ge= “dr (Xi +e) 
2 Ze 2 ee Z 
B= (Ag-f) + 4hd H=d (A t¢e)QA +f) 
2 ios 2 2 
C = (ees) + 4. d Ju="909 4)-d 
2 2 2 
D = A(Ag-e) + 2hd Ki= A> +d 
Z 2 2 2 
E = A(Ag-f) + 2hd eee eet Can ee) 
2 Z 
F = (A -e) (Ag-e) + 4dhd rn = 2(1 + Wy/Wy) 
2. 2 
a = max! J ty /2W e=d-a 
| 2.1 Ds PD 
Deseo eee at so } /2W, fea eden} 
2.45 | 
= Jed SD ary 
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Appendix III: Program to Calculate the AAX Spectrum 


Included in this section is the listing of the subprogram 
used to find the AAX spectrum. (10 tind this .spectrum, we .have to 
invert the 8x8 matrix given in table 2(b). To do this we have used 
ane cofactor method. In general, this is an inefficient method of 
solving for the inverse of an 8x8 matrix. However, we have used this 
method for three reasons. The first is that, as described in the 
comment statements of the listing, we obtain the intensity of the 
spectrum as a function of frequency when using this program. This is 
important because, for this thesis, we need the derivative of the in- 
tensity withgrespect ao kirequency. sihwsfasicasuly obtainab lemiere 
sance we havesthe antensity as: awftunct1on Or sirequencVane. 1 more usual 
methods of finding an inverse are used, one obtains an intensity val- 
ue for a certain value of frequency, the result being that one has a 
discrete set of intensity values. The fact that we obtain the in- 
tensity as a function of frequency by one inversion of the matrix 
means that computer time is saved compared to the more usual methods 
which would require that the inverse be found for each frequency val- 
me. Einally,.if,cerlain, elements, ofa matrix are zero, asampla fi; | 
cations can be made when finding the inverse by the cofactor method. 


In our case, 24 out of the 64 elements of the 8x8 matrix are zero. 
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SUBROUTINE FEED(WW,WWW,IZZZ,A54,N) 
A DESCRIPTION OF SUBROUTINES USED TO CALCULATE AA*X SPECTRUM. 


ALTHOUGH THERE ARE SEVERAL SUBROUTINES USED TO CALCULATE THE 

AA*X SPECTRUM IT IS ONLY NECESSARY TO CALL THE SUBROUTINE FEED 
AND ITS TWO ENTRY POINTS (FIRST AND SECOND) BY THE CALLING PROGRAM. 
ALL OTHER SUBROUTINES ARE CALLED BY FEED. A DESCRIPTION OF FEED 
FOLLOWS({A DESCRIPTION OF FIRST AND SECOND FOLLOW THEIR ENTRY 
STATEMENTS). 

if SHOULD BE NOTED THAT FEED IS CALLED ONCE, THEN FIRST IS CAL= 
LED 36 TIMES, THEN SECOND IS CALLED ONCE FOR EACH AA*X SPECTRUM. 


FEED; 


THE VEC- 
TORS A,ABL,DETT ARE CALCULATED. WHEN N IS NOT EQUAL TO ZERO, THEN 
THE SUM OF TERMS A(N+1)*X**N DIVIDED BY THE SUM OF TERMS DETT (M#1) 
*X**M GIVES THE INTENSITY AT FREQUENCY REAL(B) WHERE, N GOES FROM 
O TO 7; M GOES FROM 0 TO 8; AND X=-D-B, D BEING REAL AND B BEING 
IMAGINARY. TO OBTAIN THE INDIVIDUAL SPECTRA FOR THE A AND At NUC- 
LELG Nw los SET) EQUAL {TORZERO: THEN THE A AND A* SPECTA ARE OB- 
TAINED BY REPLACING THE VECTOR A IN THE ABOVE SUM EY A AND ABL 
VECTORS RESPECTIVELY. THUS BY USING THE VECTORS WITH A CALLING 
PROGRAM WHICH CALCULATES THE SUMS MENTIONED ABOVE, THE SPECTRA 
CAN BE EVALUATED. D SHOULD BE SET TO THE VALUE OF A54 UNLESS 
THERE IS BROADENING DUE TO NUCLEI NOT IN THE GROUP AA'X SUCH AS IN 
THE THESIS OF M. SMITH (1972). FOR THAT THESIS D=BJ*SJI*SJ1I*.5 
/¥W +A54 
IZZZ-THIS IS AN INTEGER. WHEN IZZZ IS NOT EQUAL TO ZERO, THE 
PRODUCT OF THE MATRIX OF THEORY AND ITS INVERSE IS PRINTED. THIS 
PRODUCT MATRIX IS NOT NORMALIZED AND THE FREQUNCY USED IN CALCU- 
LATING THE MATRIX IS ECUAL TO A54. IF 1222=0,THIS PRODUCT MATRIX 
IS NOT PRINTED. THIS IS ONLY USED AS A CHECK OF COMPUTATIONS. 
A54-THIS IS A REAL NUMBER. IT IS OF THE ORDER OF THE FREQUENCY 
FOR WHICH THE RESONANCE INTENSITY IS EQUAL TO .5 TIMES THE MAXI- 
MUM INTENSITY. 
N-THIS IS AN INTEGER. THE PUNCTION OF THIS VARIABLE IS DESCRIBED 
ABOVE. 
A DESCRIPTION OF THE PARAMETERS OF THE SUBROUTINE STATEMENT. 
WH-THIS IS THE TRANSITION RATE FOR THE A OR At NUCLEI. 
WHW-THIS IS THE TRANSITION RATE FOR THE X NUCLEUS (WW,WHW ARE IN 
1/ SECONDS) i 
IN ADDITION TO THE ABOVE VARIABLES THE FOLLOWING VARIABLES MUST 
APPEAR IN A COMMON STATEMENT OF THE CALLING PROGRAM. THE STATEMENT 
TAKES THE FORM-COMMON/ AREA/A(8), DETT (9) ,ABL (10) ,XJ1,XJ32,XJ3 
A-THIS VECTOR IS DESCRIBED ABOVE. THIS IS A COMPLEX VARIABLE OF 
DOUBLE PRECISION (COMPLEX*16). 
DETT-THIS VECTOR IS DESCRIBED ABOVE. THIS IS A DOUBLE PRECISION 
COMPLEX VARIABLE(COMPLEXK*16) 
ABL-THIS VECTOR IS DESCRIBED ABOVE. THIS IS A COMPLEX DOUBLE 
PRECISION VARIABLE(COMPLEX*16) 
XJ1,XJ2,XJ3-THESE ARE THE INDIRECT COUPLING CONSTANTS OF THE 
FIRST,SECOND, AND THIRD NEAREST NEIGHBOURS(IN HZ). 
THIRTY-SIX RECORDS MUST BE READ FROM LOGICAL UNIT 2. THE FOL- 
LOWING 36 CARDS DISPLAY THE CONTENTS OF THE 36 RECORDS, EACH RE- 
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CORD CONSISTING OF 23 LETTERS. 


BCDEFGHJKPILKLMINJIOOOO CARD1 

ICDEPGHJKPILKLMINJINMPO CARD2 

NBDEFGHILKJKPLMINJIIMPC CARD3 

MGDEBCHLIJMNIJIKLPKPINO CARDY 

PFEDBCHNMIILJKLJIPKMINO CARDS 

BCEHFGOKPILKNINJIMPJILM CARD6 

BCDEFGOJKIZTLNLMINMPPKJI CARD? 

HACFDEGKJINMPILLIMNOOOO ’ CARDS 

AHFGDEOMEI JIPLIMNJKNKIL CARDS 

KHFGADCINLJIKMLPMNIPIJO CARD 10 
AHCFDEONMIKJPILLIJKPIMN CARD11 
PEFCADGJKTILNMLNIPAKIJO CARD12 
AHGBDEFPIMIPJMNJKLIOOOC CARD 13 
IHBFAEGJLMNPJIIKMIPNKNLO CARD1I4 
LHFGDABINKJIPLAMPJIKING CARD15 
AHGBDEOPINIPKMNJKILMJLI CARD16 
KEBFDAGPJI KINJILMMPLINO . CARD17 
BCEAFGHKIPLNKINMNPJIOOOO CARD18 
BCAHFGOIPJNKIMPJILMKLIN CARD19 
LGAEBCHMIJPNIINKLPKMJIO CARD20 
MFAEBCH PNI MIJINKLPKLJIO : CARD21 
JAFGEHCLNIMPNIJNILKIKPO CARD22 
BCDAFGHJIPINKLMMPJIOOOC CARD23 
IGADBCHMLJ PAILINJIPKNKLO CARD24 
NFADBCHPMI ALJINJITPKIKLO ~ CARD25 
BCADFGOIJK NILMPLHINPKJI CARD25 
AHBCDEGINPPKIJKILMNOOOO CARD27 
ADBCEHOINM JILKPLKIJPMNI CARD28 
AHBCDEFINMPKJJKILLIOOOO CARD 29 
IECBADFKPJLKINIIJHLNPMO CARD30 
JDGBAECIPKNJIPNIKNLLAIC CARD31 
CBADFGENILIJKMPLMINOOOO CARD32 
BCHDFGOPJI KINJILMNMPKLIN CARD33 
AHCGDEONPIKIPILANJKMJLI CARD34 
AHBFDEOIMNPJKJKLIILPINN CARD35 
BCAEFGOIKJNLIMPINLMPKJI CARD36 


THE ABOVE 36 RECORDS MUST BE READ EACH TIME VECTORS ARE REQUIRED 


FOR AN AA'X SPECTRUM. AFTER THE READING OF THE 36 RECORDS, LOGICAL 


UNIT 2 IS REWOUND BY SUBROUTINE. 


COMPLEX*16 DETT,A,ABL 


COMPLEX* 16 ALL (10) , AXX, KX RET (36, 7)<g24%,DDX,X2 (40, 8) , ce 8) ,DX,DET 


$ (718) -R%_Z(7) «TT (8) 


REAL HACOS,HASIN 
COMPLEX XFAT 

COMPLEX CMPLX 

DIMENSION I1(23),III (23) 

REAL*8 AL(8,8) ,STORE (36,16) ,X4 (16) 

COMMON Z,X4 

COMMON TT 

COMMON XGAM1,XGAM2,XGAM3,XGAM4,XGAM5, XGAN6 
COMMON/AREA2/F1,XRET,STORE 
COMMON/AREA/A (8), DETT (9) ,ABL(10), XJ1,XJ2,XJ3 


DATA Iity'! ) eens SiON a UCU psa Coes Men ny of th ooUia grr pay eC ad Bs QO 
a TE, BOO IO Uy Ui toes UU ONL OY esi gO. 


DO 311 JKK=1,8 


311 ALL (JKK)=0. 


R9=-3.*WW- WWW 
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1 Aa 


308 


309 


PI=3.141593 
SPI=2.*P1T 

XJ 1=SPI¥*XJ1 
XJ 2=SPI*XJ2Z 
XJ3=SPI*XJ3 


DD=.5*SQRT (.25* (XJ2-XJ 3) ** 2+ XJ 1** 2) 


RCOS=1. 
RSIN=0. 


IF (DD.NE.0.) RSIN=. 


5*XJ1/DD 


IF (DD. NE.0.) RCOS=.25*(XJ2-XJ3) /DD 
R1= (-2. *XJ 1-XJ2-XJ3) /4.-DD 
R2= (-2. *XJ 14XJ 2+XJ3) /4.-DD 
R3= (2.*XJ1-XJ2-XJ3) /4. -DD 
R4G= (2.*XJ14X3 24+XJ3) /4.-DD 


R5=R1+2.¥*DD 
R6=R2+2.*DD 
R7=R3+2.¥*DD 
R8=R4+2.*DD 
R1=R1+A54 
R2=R2+A54 
R3=R3+A54 
R4Y=R4+A54 
R5=R5+A54 
R6=R6+A54 
R7=R7+A54 
R8=R8+A54 


IF (N.NE.0) GO TO 309 

HACOS=SQRT (.5* (RCOS+1.)) 
HA SIN=SQRT (.5*(1.-RCOS)) 
IF (RSIN.GT.0.)GO TO. 308 


HACOS=-HACOS 


IF (HACOS.EQ.HASIN)GO TC 309 
IF (HACOS.EQ.-HASIN) GO TO 309 


DUM1=HACOS+HASIN 
DUM2=HACOS~HASIN 


ALL (1) =CHPLX (HACOS/DUM2,0.) 
ALL (2) =CMPLX (-HASIN/DUM2,0.) 
ALL (3) =CMPLX (HACOS/DUHM1,0.) 
ALL (4) =CMPLX (HASIN/DUM 1,0.) 


ALL (5) =ALL (4) 
ALL (6) =ALL (3) 
ALL (7) =ALL (2) 
ALL (8) =ALL (1) 
CONTINUE 
I11=0 

DO 1 I1=1,36 
I11=11141 


READ(2,1000) (II (L1) ,L1=1, 23) 


DO 1 H1=1, 23 
42=M1-7 
IF (M1.EQ.1)M2=M1 


IF (M1.EQ.1)XRET (I11,1) =(0.,0.) 
IF (II (M1)~.EQ.III(1 ))XRET(I11,M1) =CMPLX(R9,R1) 


IF (II (M1) - EQ. III ( 
TF (il (01). £0. 112 ( 
IF (II (M1). EQ.III( 
TP (IT (M1) CEO. TI1T ( 
IF (II (M1). EQ. III ( 
IF (II (M1). EQ. III ( 
IF (II (M1). EQ. III ( 


2)) XRET (111,41) =CMPLX (R9, R2) 
3) ) XRET(I11,M1) =CMPLX (R9, R3) 
4) ) XRET(I11,M1) =CNPLX (R9, R4) 
5))X RET (111, M1) =CMPLX (R9, RS) 
6)) XRET(I11,M1) =CMPLX (R9, R6) 
7) ) XRET(1I11,M1) =CMPLX (R9,R7) 
8)) XRET (111, 1) =CMPLX (R9, R8) 
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2D 


1000 


14 


27 


IF (II (M1). EQ.III( 9)) STORE (I11,M2) =WWH*RSIN 
IF (II (M1) . EQ. III (10)) STORE(I11,M2) =WW*RCOS 
IF (II (M1). EQ. III (11) ) STORE (1I11,M2) =WWH*RCOS 
IF (II (M1) . EQ. III (12) ) STORE (111,M2) =WW*RSIN 
IF (II (M1). EQ. III (13)) STORE(I11,M2) =-WWH*RCOS 
IF (II (M1). EQ. III (14)) STORE (111,M2) =-WW*RCOS 
IF (II (M1) . EQ. III (15) ) STORE (111,M2) =-WW*RSIN 
IF (M1.NE.1)GO TO 25 

CC=AIMAG (XRET(1I11,1)) 

XYY=STORE (111, 1) 

IVUCC. 60-0.) XEET(1 19, 1) =CMELX ( XVYo 1,05) 
CONTINUE 

TERMS CT 5)GO, 10 12 

IF (11 (41). EQ. III (23) ) XRET (111,41) =(0-, 03) 
Couto 

IF (II (M1). EQ.III (23) ) STORE (111,M2) =0. 
CONTINUE 

JK=0 

RKK=2 

FORMAT (23 (A1) ) 

CALL REWIND(2) 

RETURN 

ENTRY FIRST(I11) 


HOY SHES 


I11-THIS IS AN INTEGER. WHEN FIRST IS CALLED FOR THE FIRST TIME 
FOR EACH AA'X SPECTRUM, I11=1. EACH ADDITIONAL TINE I11 IS INCRE- 
MENTED BY 1. THUS ON THE LAST CALL TO FIRST FOR EACH AA'X 
SPECTRUM 1I11=36. 
Z4=(1.,0.) 

RY=(12-02) 

JK=0 

DONS th2= 1,7 

Z (M2) =XRET (111,M2) 
CONTINUE 

DO 14 M2=1,16 

X4 (M2) =STORE(I11,M2) 
CC=AI MAG (XRET (1I11,1)) 
LE (CC 5E0 50 coer) 10 
CC=AI MAG (XRET (I11,7)) 
TE(CCLEO.0.).GOm tO | 
JK=1 

CALL ZERO(JK) 

NO wisi =178 

TT (I) =TT (I) *(-1.,0.) 
GO 10 6 

CALL NZERO (KKK) 

GO TO 6 

CALL ZERO (JK) 

GO TO 6 

CONTINUE 

DO 27 I=1,8 
¥2(111,1)-77 (1) 

Z4=TT (1) 

IF (111.GT.5)GO. TO 4 
DO°20 .1=1,8 
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DET (2.14, 1) =TT,(1) 


IF (I11.NE.8)GO TO 5 


DO, 15°1=1,8 
RX=10.,, 0.) 


IF (I.NE.1) RX=DET(1,I-1) 
DETT(I) =DET(1,1) *XRET (8,2) +RX 


DO 15 N2=2,5 
M1=N2+1 


DETT (I) =DET (M2, I) *STORE(8, M1) +DETT (I) 


DETT (9) =DET (1, 8) 
CONTINUE 

RETURN 

ENTRY SECOND 

DO, 22,1=1,8 
WG lt e204 51) 
Elie 2h (2 1) 
F1(1, 3) =X2 (3,1) 
F1(1,4) =X2 (6,1) 
FI.(1,6) =x%2 (8,1) 
FILO, je %2 (5,2) 
Fal (lp 8X 23(7,, 1) 
F1(2,2)=X2 (8,1) 
F1 (2,3) =X2 (9,1) 
F1.(2,5) =x2.(10, 1) 
F1(2,7) =X2 (11,1) 
F1 (2,8) =X2 (12,1) 
F1 (3,3) =X2 (13,1) 
F1(3,4)=X2 (14,1) 
F1(3, 5) =X2 (15,1) 
F1 (3,6) =X2 (16,1) 
F1(3, 8) =X2 (17,1) 
F1(4,4)=X2 (18,1) 
F1(4, 5) =X2 (19,1) 
F1(4,6) =X2 (20,1) 
F1(4, 7) =X2 (21,1) 
F1(2,4) =X2 (22,1) 
F1(5,5) =X2 (23,1) 
F1(5,6)=X2 (24,1) 
F1(5, 7) =X2 (25,1) 
F1 (5,8) =X2 (26,1) 
F1(6, 6) =X2 (27,1) 
F1 (6,7) =X2 (28,1) 
F147; 7) =X2 (29,1) 
F1 (7,8) =X2 (30,1) 
F1(6, 8) =X2 (31,1) 
F1 (8,8) =X2 (32,1) 
F1( 1,5) =X2 (33,1) 
F1 (2,6) =X2 (34,1) 
F1(3, 7) =X2 (35,1) 
F1 (4,8) =X2 (36,1) 
DO 13 N2=1,8 

DO 13 N1=1,8 


F1(N1,N2) =F1(N2,N1) 
TEMIENE 1) GO TO 21 


CALL CHECK (1222) 
D1=1.-RSIN 
D2=1.+RSIN 
D3=RCOS 

AL (1, 1) =D1 
AL(1,2)=-D1 
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300 


301 


304 


302 
21 


DO 300 J1=3,6 

AL (1,31) =-D3 

AL (1,7) =D1 
AL(1,8)=-D1 

AL (3, 1) =-D3 

AL (3,2) =D3 

pO 307 J1=3,6 

AL (3,91) =D2 

AL (3,7) =-D3 

AL (3,8) =D3 

DO 304 J2=4,6 

DO’ 308 J1=1,8 
AL(J2, 01) =AL (3,31) 
DO 302 JI=1,8 
D4=AL(1,J1) 

AL (7,31) =D4 

AL (8,31) =-D4 
AL(2,J1)=-D4 
CONTINUE 

CONTINUE 
A(I)=(0.,0.) 
RBU(T) = (04, 0%) 

DO 303 J1=1,8 
AXX=ALL (J1) 

DO 303 J2=1,8 

IF (N-NE.0)GO TO 305 
ABL (I) =AL(J1,32) *F1(J1,32) *AXX+ABL (I) 
A(I)=AL(J1,32) *F1(J31,J 2) +A (I) 
CONTINUE 

CONTINUE 

IF (N.EQ.0) GO TO 306 
RETURN 

CONTINUE 

IF (ABS (HACOS) -NE.HASIN)GO TO 313 
pO 314 JKK=1,8 

ABL (JKK) =A (JKK) /2. 
CONTINUE 

DO 307 I=1,8 
A(I)=A (I) -ABL (I) 
RETURN 

END 


SUBROUTINE NZERO (KKK) 


REAL*8 Heyl eJe Ke le Me Ne Pe Qe ReSeT,UeV,W,X 
REAL*8 XGAN1,XGAM2,XGAH3,XGAN4,XGAN5,XGAM6,A1,A2,A3,A4,A5,A6,A7, 


$A8,A9,A10,411,A12,413,A14,A24,A16,A17,A18,A19,A20,A21,4A22,A23,A 


COMPLEX* 16° Y,8,C,D,&,F,G, XX (8) -TTs ZZ (8) <C1,.C2 

COMNON YB TCD EECE SG Held, Kiel, My Ny Pee ReS, fy U s.V, REX 
COMMON 2Z 

COMMON XGAM1,XGAM2,XGAM3,XGAM4,XGAM5, XGAM6 
COMMON/AREA1/C1(10, 8) ,C2(10,8) ,XX,TT 

ary 

XGAM5=N *R- P*Q 

XGAM6=N*T-E*S 

A1=H¥*U-K*A 
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A2=R¥*S-Q4T 

A3=M*R-L*¥T 

A5=L*S-4 *Q 

A6G=N*M-K*S 

AT=M¥A-J 4U 

AS=E*M-K*T 

AQ=S#W~VAT 

A10=-A2 

A11=Q*W-V*R 

A12=W*L-U*E 

A13=W*M-T*U 

A14=a9 

A16=W4*N-V*P 

A17=L*P-K* E 

A18=I¥M-L*J 

A19=I1*U-L*A 

A20=L*N-K*Q 

A23=K*J-M*H 

A24=I*K-L*H 

TTsA{0 405) 

DOo5ouli=t-c 

IX (11) = (055.05) 

Docdeti= 1710 

DO, 4 JJ=1,.8 

C2 (LE JJ (0 400-) 

CIA il id) = 10,05) 

C1(1,1)=F 

C1(2, 1).=D 

C1(3,1)=E 

C1(4, 1) =c 

C1(5,1).=G 

C1(6, 1) =B 

DOs? 1-156 

CUCL c= Gln Oe) 

DO 9 II=1,8 

XX (25) = (05,0) 
CALIGGALC (2, 3,2, 10} 

pO 10 I1=1,8 

C110 ,11)=C2 (10, 21) ey 

DOV ities 

XX (11) = (0..,0-) 
XGAM1=-K¥*A 2-L* XGAM64M* XGAM5 
XGAM2=U*XGAN6+V*A8-W*A6 
XGAM3=-U*A 24+ V*¥A3¢ WAS 
TT=A2*K¥*(A1*A2 #A*(P*AS+#N*A3) +H¥* (M¥A11-AO*L)) +XGAMS* (M*A*(2.*L* 
$ XGAM6-R*A6) +A3* (J*®N*U +H*V¥M) +0*S *U*A1740%A 7*A 8) + XGANG*L* (A*(A17*S 
$-T*A20) +I*U*¥XGAMNG41* (K *A14-M4*A16) ) +AZ¥N* (OX ( -A13*H) +A12*H*S) 
TT=TT+ XGAM1* (J* (K¥A11-L*A16 ) ~-XGAM2*1) ~1*L* XGAM2* KGAM6+ 
SXGAM3*S*L* P#H+M*P¥Q*H* (-W*A5+ UA 2) 
TT=TT+A3*(-H*L4*N¥*V*T) 

XX0(1)=1T 

TT=AG6* (S¥*A 14N*A7+V*A 23) 
CALIOCAUC Glee ped) 

TT=A9*(T*A 14#P¥*A7+9*A23) 

CRI eChIC(3,2,.2,2) 

TT=-A5* (S*A194Q*A7-V*A 18) 

CALI CALC (1, 4,2, 3) 

TT=A3* (T¥A19+R*A7-W*A18) 

CALL CALC (3, 8,2, 4) 
TT=A20* (OX A1-N*A19+V*A 24) 
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CALL CALC(1,5,2,5) 

TT=A2* (A*A 2-1 *A14+J*A11) 
CALDZCALCHG 4 ,2,6) 
TT=A17* (R¥A1+W*A24-P¥A 19) 
CALIVCALC(3,5,2,7) 
TT=XGAM6* (A*¥XGAM6+H*A14-J*A 16) 
CALL CALC(6,2, 2,8) 
TT=XGAM5* (XGAMS*A+H*A1 1-1*A16) 
CALL CALC(5,6,2,9) . 
DO 5 II=1,9 

DO 5 JJ=1,8 

©€4 411533) =C2 (11,34) 
TT=A19*L 
CALISGALC(3,7,4,2) 
TT=A1*K 

GALIICALC (1,754,.1) 

TT= (J *V-S*A) *S 
EALIGEALC( 1,6; 473) 

TT= (J *W-T*A) *T 
CALLACALC(2,6,4,8) 

TT= (V*I-Q*A) *Q 

CALL CALC (5,6,4,5) 

TT= (I*W-R*A) *R 

CALL, CALC (7,6,4,6) 

TT= (H¥V-N*A) *N 

CALIVGALC (5,074 ,7) 
TT=-A7*M 

CALL CALC (2,3,4,8) 

TT= (W*H-P*A) *P 

EAT CALC(9;2,4,9) 

DOSE TI=17e 

C1 (511) =C 25, 22) 
TT=(12,02) 

CALI CALC (5,10,6,1) 

DOw 2et1=eA, 6 

ZZ (IL) =XX (II) * (-1.,0.) 
RETURN 

END 


SUBROUTINE ZERO (JK) 


REAL*SuH DJ, Kecle Me Neh eGeReSetelevews x 

GOMNBLESE AG C1, C2, X17 k2,%3, 04 kop kOpkk (Sette lo Olbs ke Br Cs DE, E.G 
REAL*8 XGAM1,XGAM2,XGAM3,XGAM4Y,XGAM5, XGAM6,A1,A2,A3,A4,A5,A6,A7, 
$A8,A9, A10, A11,A12,A13,A14,A15,A16,817,A18,A19,A20,A21,422,A23,A 28 
REAL*8 A25,A26,A27,A28 ,A29,A30,A31,A32,XGAM7 

COMMON A,B Ce, D,EvEsG;Hrlsed okeleMeNe PeQe Be Sete Ve Wek 

COMMON ZZ 

COMMON XGAM1,XGAM2, XGAM3,XGAM4,XGAN5S, XGAN6 
COMMON/AREA1/C1(10, 8) ,C2(10,8) ,XX,1T 

XGAM5S=P*C-N*R 

A1=T*N-P¥*S 

AS=C*¥T-R*S 

A15=L*T-M*R 

A16=N*P-K*T 

A17=L*S-N¥*Q 
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A18=M*N-K*S 
A22=H*S-J¥*N 
A24=P*J-H*T 
A26=I*T-R*J 
A28=J*K-H*Ms 
A30=I1*H4-J¥*L 
A32=1*S-J*¢ 


IF (JK.EQ.0.)GO TO 1 


U=J 


(Hn 


4 


& Kh we < 


=-A5 
A6=-A15 
A7=A17 
A8=-A18 
AS=A16 
A10=A1 
A11=-A26 
A12=-A32 
A13=A24 
A14=-A22 
A23=A13 
A25=A26 
A29=A28 
A31=A30 
CONTINUE 
A2=L*P-K*R 
A3=K*Q-N*L 
A19=R*H-P¥*I 
A21=Q*H-I*N 
A27=H*L-1I*K 


IF (JK2EQ21)GO TO 2 


A4=R*¥WH-Q*¥X 
A6=V*R-L*¥X 
A7=L*W-V¥*Q 
AS=K*W-V¥*N 
AQ=V*¥P-K*X 
A10=X*N-P*W 
A11=R*U-I* xX 
A12=U*Q-1*W 
A13=U*P-H* xX 
A14=U*N-H*W 
A23=U*P-H*X 
A25=I1*X-U*K 
A29=U*K-H*V 
A31=1*V-U*L 
CONTINUE 
A20=-A10 


XGAMN1=-U*A Z-H*A6+I¥*¥A9 
XGAN2=H*A7 -I*A8+U *A3 
XGAM3=-Q¥*A 13-I1*A104R¥*A 14 
XGAM6=—-J*A2+H¥*¥A 1541 *A16 
XGAM7=H*A17+I*A18+J*A3 


TT=(0.,0.) 
DO 58 I1=1,8 


XX (11):= (0/702) 

DO. 20 siT=1 710 

DO 20 JJ=1,8 

Gtld 3455.3) = 10.0.) 
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CQUTESII)= (02,0) 
C1(1,1) =c 

C1(2,1)=D 

C1(3, 1) =A 

C1(4,1)=B 

C1(5,1)=E 

C1(6,1)=F 

C4172 1;=c 

DO"Z1 FI 1=1 710 
G441ts2)=(1.,0-) 

CALI CALC (3,4,2,1) 
FALISCALC( 653, 262) 
EALLIGALC (3557273) 
CALLIGALC(6, 4, 2,4) 
CALIGECAEC (5,4, 2,5) 
CALEICALC(6, 5, 2, 6) 

CALI, CABG(3;272;,9}) 

GALL CALC (6,252, 10) ; 
X1T=L*¥A1* (W*A2+X*A3-V*X GAM5S) +A5*K* (K¥AG-N*¥AG-P¥AT) +M*#XGANS* (-R¥A8B 
$-Q* A9-L*A10) 
X2=XGAM3*XGAM5*J+H*A5* (H*¥AG-N*A114P¥*A12) +I*¥A1* (-I1*¥A10-Q*A134R*A 14) 
X3=-XGAM1* XGAM6 
X4=-XGAM2*XGAH7 

DO 52 I1=1,8 

XX (11) =(0-/-0:) 

XX (1) =X1*A+X2*B+XZ3*EtX 4*P 
XX (2) =X 14#X 24X34X4 
X1=-A20*A1 

X2=-A18*A8 

X3=A16*A9 

X4=-A22¥*A14 

X5=A23*A24 

X6=A28*A29 

DO 26 JJ=1,8 

C1 (10,3J) =X1*C2 (1,55) #X2*C2(2,5I) #X3*C2(3,5I) #X4*C2(4, TI) #X5*C2 
$ (5,J3I) +X6*C2(6,JdJ) 

T= 6149202) 

CA14 GALG(2,10 23,7) 
X1=-AY#*A5 

X2=A17*A7 

X3=-A6*A15 

DO 27 JJ=1,8 

C1(10,J3J) =X1*C2(1,JJ) #X2*C2 (2,53) #X3*C2 (3, JJ) 
TT=(1.,0.) 

CALE GALC(1¢10 ¢3¢7) 

DO 7 II=9, 10 

DO 7 JJ=1,8 

C1 (21 JJ) =C2 (i100) 

DO 8 II=4,6 

DO 8 JJ=1,8 

CALMED 5)=C2 (11,33) 
TT=-A32*A12 

CALL CALC (174, 3,8) 
TT=A25*A 26 

CALL CALC (1,5,3,9) 
TT=A30*A31 

CALL CALC(1,6,3,10) 

DO, 9 Tt=15:8 

C1 (3,211) =C2 (9,11) 

X1=M4V 
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11 


12 C1 (4,11) =X 1*C2 (3, IL) +X 2*C2 (5, ITI) +X 3*C2(2, 11) #X4*C2 (4,11) +X5* 
$C2 (1,11) +C2(6,11) *X6 


13 
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23 


X2=S*W 

X3=T#*X 

X4=U¥I 

DO 10: IT=1,8 

Can(W, iT )=C2:(10 ,11) *X14+ C2 (8,11) *x2 
TT=4~1.,0-) 

CALIRCALC (N29. 5 ei) 

DO 11 II=1,8 

C1(5, II) =C1(9, 11) *X3+C 1(10,I1) *X4 
TT=(-1.,0.) 

CALUCCALC (375,155 7) 

TF (IK £0.0).G0 10.22 

X1=A2*A2 

X2=A19*A19 

X3=A3*A3 

X4=A21*A21 

X5=XGAN5*XGAM5 

X6=A27*A27 

DO 12 II=1,8 


TT= (1. ,0<) 

CHUL MCA UC (437178 £7} 

DO 13 II=1,8 

ZTE) XX (TT) 

X1=N*N 

X2=E*P 

X3=K*K 

DO 14 II=1,8 

C1(5, 11) =X 14C2 (4,11) #X2*C2 (5,11) +X3*C2 (6,11) 
CALTOCAUG(5 9, 407) 

X1=R*R 

X¥2=Q*Q 

X3=L*L 

DO 44 II=1,8 

Ci(3501) = (02,04) 

C1(3,1)=A 

CVUIE2)E (Us 20".) 

DO 15 II=1,8 

C1 (5,11) =C2 (9,11) *#X 1#X2*C2 (8, 11) #C2(10, II) *X3 
CALL CALC (5,3,4,6) 

CALL CALC(9,1,3,4) 

DO 38 II=1,8 

€1(10, T1j=]=C2(5,11) 

€1(6,1)=F 

CL CeR2=( 0 sh0.) 

€1.(6,53) = (02703) 

CAI, CAUC@(E9 10 433110) 

DO 4O II=1,8 

Ct (CAG AIDS C210, 1) 

X1=H*H 

X2=1*1 

DO 16 II=1,8 

C1(5, II) =(-1.,0-) *C2 (4,11) #C1(2,IL) *X1#C1(1,11) *X2 
CRLET, CALE (5 21056, 5) 

po cl Sitates 

C4 CU,TTySC2(5R1I) +C2(7 , IT) +C2 (6,11) 
DO (2341 r=1 76 

XX (II) =ZZ (II) 

TT=(-1.,0.) 
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CALES CALCACGE7 47,5) 
CONTINUE 

DO 25 II=1,8 

ZZ (1I}=XX (13) #(- 1.70.) 
RETURN 

END 


SUBROUTINE CALC(L,K,M, KK) 


COMPLEX*16 C1,C2,XX,TT 
COMMON/AREA1/C1 (10,8) ,C2(10,8) ,XX (8) ,TT 
DO 4 I=1;8 

C2 (KK, T3210 2502) 

KR=N+2 

DO 2 N=2,KR 

DO “o1=A;8 

DO 3 J=1,8 

IF(1I.GE.N) GO TO 2 

IF (ItJ.NE.N)GO TO 3 

C2 (KK,N-1) =C2 (KK, N- 1) #C1(L,d) *C1(K,I) 
GOP TON 

CONTINUE 

CONTINUE 

XX (N-1) =C2 (KK, N-1) ¥TT+ XX (N-1) 

RETURN 

END 


SUBROUTINE CHECK (TI) 


COMPLEX*16 DX 
COMELEX*8 DDX 
REAL*8 STORE 
COMPLEX*16 F1,XRET, F2(8,8) 
COMMON/AREA2/F1(8,8) ,XRET (36,7) , STORE (36,16) 
IF (I-EQ.0) GO TO 6 
DO 12 N1=1,8 

DO 12 N2=1,8 
F2(N1,N2)=(0., 0.) 
F2 (1,1) =XRET (8,2) 
F2(1,2) =STORE (1,4) 
F2 (1,3) =STCRE (1,10) 
F2 (1,6) =STORE (1,8) 
F2 (1,7) =STCRE (1,3) 
F2 (2, 4) =STORE (1,1) 
F2 (2,5) =STCRE (1,2) 
F2(2,8) =F2 (1,7) 

F2 (3,4) =F2 (1,2) 
F2(3, 5) =STORE (1,5) 
F2 (3,8) =F2 (2,5) 

F2 (4,6) =F2 (3,5) 

F2 (4,7) =F2 (1,6) 

F2 (5,6) =F2 (1, 2) 
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F2 (5,7) =F2 (1,3) 

F2 (6,8) =F2 (2,4) 

¥2 (7,8) =F2 (1, 2) 

F2 (8,8) =XRET (1,7) 

DO 16 N2=1,6 

N3=N2+1 

F2(N3,N3) =XRET(7,N2) 
DO 13 N2=1,8 

DO 13 N1=1,8 

F2 (N1,N2)=F2 (N2,N1) 
DO 14 N1=1,8 

DO 14 N2=1,8 

DXK= 105, 0-2) 

DO 15 N3=1,8 
DX=DX4+F2(N1,N3) *F1(N3, N2) 
DDX=DX 

CX=AIMAG (DDX) 
EX=REAL (DDX) 
WRITE (6, 2002) CX,EX 
FORMAT(* *,£14.7,3X,E14. 7) 
RETURN 

END 
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Appendix IV The Effect of More Distant Nuclei in the Ruderman- 


Kittel Model 


We use the approximate lineshape of an XA configuration for 
the determination of the effect of more distant nuclei. If Jay << Wy 
the XA, spectrums Lorentzian: and) 1s equal to), (43)) ofjsection i ex- 
eeptyrnac 2Wy is replaced by 

Z 
ae + ieee poy (1) 
(D.G. Hughes, private communication). To take into account the effect 
of the more distant nuclei we assume that the apparent increase in the 
relaxation rate of the observed nucleus due to each of the more dis- 
tant unlike magnetic nuclei is equal to 7(3,.1) /40W, and that the 


amcteases dd linearly. [his tis equivalent to assuming tual, thesmore 


distant nuclei are independent of each other. Then the apparent in- 


crease in the relaxation rate of a Sn!!7 nucleus is 
(7/40W aca 0 
Xo aWGin AP nbtg (2) 
36 te Ds 
= (7/40W a 2 oF (k IR. 1)p Fat TO5Sx<10 6 Sani eotey Seyt lm 0 amg ee" 
X DBA ae aes un Hf erin) 


Here the sun) S oven al lmnuclear, sites-outsidcmne sam Inncres nels. 
1 


Also w is the magnitude of the Ruderman-Kittel interaction for near- 
est neighbours, F(x) is given in equation (34) and p- is the proba- 
un 


117 nucleus) being found 


bility of a magnetic nucleus (other than a Sn 
at a lattice site. The result is that the apparent increase in the 
relaxation rate is 


2 
0.305457 / 2H, (3) 


We have not taken into account the effect of like nuclei nor the mu- 
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tual interactions between the more distant nuclei in the calculations 
Jeading to (4). Because (3) “isan approximate expression for the 
effect of more distant nuclei, we multiply this. expression by & 
and obtain for the apparent increase in Wy 

E{.30545,°/2N,) (4) 
Here € is a fitting parameter and its magnitude is found by mini- 


muzine, (2) of sectron WV. 
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